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ABSTRACT 

Radiation energy transfer and ionization rates in the mixture of 
atoms, ions and electrons produced by a normal shock wave propagating in 
cold argon gas are examined. It is shown the gas may be regarded as a 
combination of two gases, one composed only of electrons and one com- 
posed of both atoms and ions. Temperature differences between the elec- 
tron gas and the atom- ion gas significantly affect radiation and ioniza- 
tion rates in the shock wave. 

The analysis uses the ionization rates due to both atom-atom 
collisions and electron-atom collisions. In addition, photoionization, 
which is responsible for precursor ionization, is included. 

Although argon is not a grey gas, it is shown that the radiation 
variables may be reduced to the forms they would have if the gas were 
grey, but with different source functions and attenuation properties. 
Radiation generated by photo ionization and radiative recombination 
involving excited atomic states escapes ahead of the shock wave and 
results only in cooling the hot gas. Radiation due to photoionization 
and recombination involving the ground state is trapped by the cold 
particles ahead of the shock and causes precursor ionization. Photo- 
excitation processes (line radiation) are omitted from consideration. 

The integral form of the radiation variables causes the equations 


to be integro-differential. Iterative techniques for solving the 
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equations are not practical because they converge too slowly. Instead, 
a combined perturbation-iteration method is used. 

Hie effect of viscosity and heat conduction is to d is continuously 
change the temperature and density of the gas. Preceding the discon- 
tinuity is the precursor region caused by ground state photoionization. 
After the discontinuity the electron gas is colder than the atom- ion 
gas . Subsequently the electron gas temperature increases until it 
becomes equal to the temperature of the atom- ion gas. Finally, after 
chemical equilibrium has been reached, the gas radlatively cools to its 
final stage. The ionization rate and precursor ionization are found to 
be greatly affected by the cool layer of the electron gas following the 
discontinuity. 

Radiation energy transfer has little effect on the ionization rate 
until equilibrium ionization is approached. The net collisional rates 
then decrease and the radiative rates be com relatively more important. 
The net effect is to decrease the ionization rate and delay chemical 
equilibrium. 
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CHAPTER I 

BACJGQROUHD AND LITERATURE SURVEY 

In the lnet few yeers there has been e greet deel of Interest 
concerning shock waves propagating with sufficient speed that the shocked 
gas becomes ionised and emits radiation. Shock tube research has shown 
how elastic and inelastic collisions between electrons, ions, and atoms 
of the gas influence the ionization process. The role of radiation la 
the ionization process, however, has never been determined. 

In this thesis the effect of radiation energy transfer on the 
ionization process in chock waves propagating in argon is investigated. 

It is an especially difficult problem because the radiation energy trans- 
fer, the .kinetics of ionization and the energy exchange between electrons, 
atoms and ions in the gas must be considered simultaneously. The fre- 
quency dependence of the radiative emission and absorption processes 
leads to additional difficulties in calculating the radiation energy 
transfer. 

Although the analysis in this thesis is performed for shock 
waves in argon, its concepts and techniques are generally applicable 
to chemically reacting, radiating gases. 

All shock waves analyzed in this thesis are assumed to be normal 
and time independent. In the Bulerian frame of reference the three 
appropriate fluid conservation equations are: 
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The Conservation of Miss Bguetion. 

| (p “) * 0 (!•*) 

In this aquation p is the aass density of the gas, u Is the velocity 
of an eleaental gaseous volume and x is a coordinate parallel to the 
direction of shock wave propagation. The origin of x is fixed to the 
shock wave. Values of x increase in a direction opposite to the direc- 
tion of shock wave propagation as shown on Figure 1.1. 

The Conservation of Moment tan Equation, 

[pu 2 + p - u] « 0 (1.2) 

Here p is the hydrostatic pressure in the gas and a is the compres- 
sive stress produced in the gas because of its viscosity. 

The Energy Conservation Equation. 

~ [pu (h + £u") - ou + q } .0 (1.5) 

The quantity h * (e + p)/p stands for the enthalpy per unit mss of 
gas, e is the internal energy of the, gas per unit volume and q is the 
net beat flux in the direction of increasing x. 

These equations may be integrated once to give 


pu «= M 

(i.*0 

Mu + p - o s» P 

(1.5) 

M (h + ju c ) - ou + q » E 

(1.6) 


where M, P and E are integration constants. These are the equations 
which are applied to the analysis of shock wave structure in the simplest 




Figure 1.1. Coordinate System. 
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cases. It Is necessary to modify then vhen the shock wave occurs in a 
partially ionized gas. This Modification is discussed in Cfcapt** II* 

The character of the flow in a shock wave is determined by the 
quantities h, p, j and q. The way In which these quantities are In- 
fluenced by viscosity, chemical reactions, radiation and ionization is 
discussed in the following sections. 


Effect of Viscosity 

The viscosity of a gas is responsible for the compression and 
resultant heating which occur in shock waves. At low speeds where the 
gas is not hot enough to produce changes of state in the gas or emit 
radiation, it is the only process which occurs. The pressure, enthalpy, 
compressive stress and heat conduction are then described by the follow- 
ing familiar relations. 


pRT 

( 1 . 7 ) 

C p T 

(1.8) 

du 

dx 

( 1 . 9 ) 

-k « 
c dx 

(1.10) 


Far upstream and far downstream the compressive stress and heat 
flux are both zero and the initial and final states are related by 

p u = p + u + * M (1.11) 

Mu + p = Mu + p (1.12) 

— • «r + 


h-+ i U_ 


V A 2 

h + Au 
+ 2 + 


( 1 . 13 ) 
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Subscripts ♦ sod - in tbs so equations denote values at + w and - m 
on the x coordinate. 

Mathematically speaking the shock extends from - » to + ». In 
practice the thickness of the viscous shock is often negligibly small. 

It may then be thought of as producing an instantaneous change in the 
state of the gas with the initial and final states related by equations 
(1.11) through ( 1 . 15 ) . In mathematical terms the change of state produced 
by the shock wave may be described as a step function. 

Effect of Finite Chemical Rates 

The final state reached by the gas in a shock wave can be chem- 

ically changed compared to its initial state. Description of chemical 

changes in shock waves can be very complicated. For example, Sherman 
£ 

[1] has studied a radiating shock wave in air involving thirty-four 
reactions. The shock waves studied in this thesis are assumed to propa- 
gate in argon* . The choice of argon is advantageous because the only 
chemical process which can occur in argon is ionization. Consequently 
the chemistry is relatively uncomplicated. Furthermore, because of its 
wide use in shock tube experiments, the properties of argon are v#ell 
known. 

For ionized argon gas the enthalpy and pressure are given by 
the following formulas [2], 

h = | (l + a) pRT + a R T lon (1.14) 

P « (1 + a) pRT ( 1 . 15 ) 

Numbers in brackets refer to the Bibliography. 
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In these equations o is the degree of ionisation in the gas 


a - 


“a ♦ n l 


( 1 . 16 ) 


where is the number density of ions and n ft is the number density 
of atoms. T l0Q is the ionization temperature defined by the require* 

* • 

ment that kT ion be the ionization potential. For argon T ion » 1.82(l(r ) K. 

If the ionization is rapid, the gas may be thought of as instan- 
taneously adjusting to a degree of ionization such that the rate of 
ionization and the rate of recombination cancel each other. The gas 
is then said to be in chemical equilibrium. If chemical equilibrium is 
maintained in the gas, the degree of ionization may be calculated from 
the Saha equation [2] which for argon may be written 


2 

Q 

1 - a 


ion 


(-S-)' 

\ x ion / 


3/2 


- T 


ion / 


(1.17) 

where P^ on Is 8 constant which has the value 149.3 gm/cm^. 

However, because of the finite rate at which ionization processes 
proceed, the changes which would take place in a shock wave if chemical 
equilibrium were maintained in the gas, may be more rapid than can 
actually occur. The degree of ionization in this case must be calculated 
from equations which describe the chemical kinetics of the ionization 
process. Assuming no diffusion of electrons, conservation of electrons 
requires 


~ (n u) m h 
dx v e ' e 


+ n 


'coll rad 

where n g is the number density of electrons, & 


e 


coll 

production electrons by collisions! processes and & 


( 1 . 18 ) 

is the net rate of 
'rad 


is the net 
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rate of production of electrons by radiative processes. If the nunber 
density of electrons and ions are equal, the degree of Ionisation may be 
written 


n. 


“a + “t 


“a + “i 


n e 

7T 


where is the atomic mass. Equation (1.18) may then be transformed 

to 


M 


oa 

dx 


m (ft 

M ' l 


a ' e 


+ A 


) 


'coll rad 

Collisional ionization occurs by the following reactions; 


(1.19) 


A(p) X ~ A + + X + e p = 1,2, ... (1.20) 

A(p) + X t A(q) + X q> p, p » 1,2, ... (1.21) 

A represents an argon atom, e an electron and X stands for a typical 
reaction partner which could be A, e, or a contaminant species. Bound 
electron states are indicated by the index p or q with p = 1 corresponding 
to the ground state. 

Radiative ionization is described by a similar set of reaction 
equations. 


A(p) + hv A + + e 

p = 1,2,5, ... 

(1.22) 

A(p) + hv «- A(q) 

q> p, p * 1 , 2 , 5 , ... 

(1.23) 


Continuum radiation is emitted and absorbed by reactions (1.22) and 
line radiation by ( 1 . 25 ) • 

Other reactions are possible but the research to be reviewed 
presently shows them to be of secondary importance for the calculations 
performed in this thesis. 


L 
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Detailed expreaaions for the coillslonal and radiative ionization 
ratea are developed in Chapter XI. 

Biffect of Radiation 

Radiation of the gas in a shock wave can affect shock wave struc- 
ture by contributing to the ionization process, as Just discussed, and 
by its contribution to the heat flux q. The effect of radiation on 
pressure and energy density is completely negligible for the calculations 
performed in this study [5]. 

The equation of radiation energy transfer appropriate to a one- 
dimensional geometry is 

u Sr <s v) ■ - — i — ll(>, u, v) + 3(x, v ) ] (1.SA) 

1 (x,v) 

The radiation intensity I(x, u, v) is defined by the requirement that 
I(x, u, v) dA dto dv is the radiation energy which passes through area 
dA of a surface at x whose normal is parallel to the x axis, within 
solid angle do) of a direction making an angle cos’^p with the x axis, 
in the frequency interval v, v + dv, per unit time. The geometry approp- 
riate to this definition is pictured in Figure 1 .2. 

The source function S(x, v) is the radiation intensity at position 
x, angle cos”^p and frequency v due to radiation emitted per unit mass 
of gas at the location of dA. Since emitted radiation is independent of 
direction, the source function has no dependence on p. 

The penetration length i(x, v) is the length required for the 
radiation intensity to decrease by a factor l/e in a gas of constant 
properties if there is no emission of radiation (as can be seen directly 
from equation (1.24)). The penetration length Is a function of x 
through its dependence on temperature, density, etc., and it is a faction 


h 


SHOCK 


WAVE 



Figure 1.2. Radiation Transfer Geometry. 




10 


of frequency through the way in which radiation is created and destroyed 
by the atoms, ions and electrons in the gas. Detailed formulas for the 
penetration length are developed in Chapter II. 

At the upstream boundary x ■ x_ the gas is assumed so cold that 
no radiation is emitted in the positive x direction. At the downstream 
boundary, the gas is assumed to have constant properties and the inten- 
sity is given by the Planck function B v (T ) ■ 2hvVc^(e* 1V// ^-l)"^ 
evaluated at the upstream temperature T + . 

I(x^, n, v) *o for 0 a u z l 

( 1 . 25 ) 

I(x + , u, v) - B v (T + ) 


Replacing x by the optical length 
x 

t(x, v) = J 


dx 1 


( 1 . 26 ) 


,V) 


in terms of which equation (1.2k) is linear, the solution of equation 
( 1 . 2 k) appropriate to the boundary conditions ( 1 . 25 ) is seen to be 

C t ( x > v ) 

J s(x',v) exp lr [t(x, v) . r(x', V) )/p)dT(x', v)/n 

for 0 s u s 1 


t(x.,v) 


I(x,u, V) « * 


B V (T + ) e 


-(t(x,v) - t(x + ,v)]/u 


( 1 . 2 ?) 




t(x + ,v) 

- J s(x,v) exp {-1 t(x, v) - f(x», , v)/u 

t(x,v) for -1 * u i 0 
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The contribution to the heat flux from radiation in the frequency 
range v, v + civ ia 


(^(x, v) - 2 kJ u I(x, **, V) dn 


( 1 . 28 ) 


Substitution of (1.27) into this expression gives 

t(x + ,v) 

q R (x, v) - 2« I J S(x’v)86n(T(x,v) - t(x’,v) ]E £ ( |t(x,v) - t(x* , V) | )cLt(x' # v) 


- e 5 [t(x + ,v) - T(x,v)) B v (T + ) 


where 


A 

Vt) * J u n * 2 e‘ t/4 du 


(1.29) 


(1.50) 


sgn (t-t* ) 


{ 


+1 if t' < t 


-1 if t« > t 


(1.51) 


The location of the downstream boundary x will be taken to be 
- oo and that of the upstream boundary x will be taken to be + •». Then 

■f 

equation (1.29) becomes 
+00 

( 

q R (x, v) • 2 % J S(x* ,v)egn[T(x,v) - t(xJ v ) iBgl |t(x,v) - r(x' ,v) | ]dT(x» ,v) 

-00 

(1.52) 

The total contribution to the heat flux due to radiation of all fre- 
quencies is found from this by an integration over frequency. 

00 

q R U) a J ( l R ( x > v ) dv ( 1.55) 

o 

The form of the radiative heat flux q_ considerably complicates 

K 

the problem because of the frequency integration required in equation 
(1.55) and because the integrand in equation (1.52) depends cm the 
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tomparaturs, density, etc. of the gti which art not known until the prob- 
l«m la solved. In alaoet all of the work on radiation shook waves tbs 
problea has boon siaplified by assuaptions involving ana or both com- 
plicating factors. 


kffects of Ionisation 

Ionization in shook waves Is responsible for some phenomena which 
are mentioned here and discussed in greater detail later. 

A partially ionized gas is a mixture of atoms, ions, and elec- 
trons. In such a mixture it is possible for processes to occur in 
which one species is preferentially affected. 

For example, electrons will lose energy by creating ions through 
the reaction 

e + A(p) A + + e + e. 

It happens that the ways in which electrons can gain energy are much less 
efficient than the way they lose energy by creating ions. The net result 
is an energy drain from the electrons. The electrons can be thought of 
as cooler than other species when ionization is taking place. 

The thermal conductivity of electrons is much larger than for 
other species. Consequently a preferential heat conduction in the 
electron gas can be expected. 

Transport properties are greatly affected by ionization. Vis- 
cosity and heat conduction are much different for a slightly ionized 
gas than one which is fully ionized Ik]. 

Collision processes by which energy and momentum exchange take 
place also depend on the degree of ionization in the gas. If the 
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degree of Ionization is snail, energy and momentum exchange takes place 
through binary collisions. But when the degree of ionization is large, 
energy and momentum may be exchanged by interactions involving many 
particles. 

Lastly, the steep concentration gradients in a shock wave will 
promote diffusion of electrons relative to Iona. The resulting charge 
separation will cause an electric field which can affect the flow. 

Literature Survey 

The survey presented in this section does not attempt to give an 
exhaustive account of all the research conducted on radiating, chemically 
reacting shock waves. Rather, its purpose is to illustrate the experi- 
mental and theoretical knowledge which has been gained by such research. 

Radiation in Shock Waves 

All of the work involving radiation energy transfer in 6nock 
waves, with the exception of the work of Clarke and Ferrari [5], has 
assumed the gas to be grey and the radiation to be emitted in local 
thermodynamic equilibrium. 

A gas is grey when its penetration length may be replaced by an 
average penetration length which is independent of frequency. 

i(x, v) £ avg (*) (1*34) 

local thermodynamic equilibrium means that the source function at 
any point in the gae is given by the Planck function evaluated at the 
local gas temperature. 

S(x, v) - B y (T) 


( 1 . 35 ) 


>1! 


li» 






With thtse simplify lag assumptions equation (1. 53) for the radiative 
flux becomes 

f + m * 

q R “ 9 J y T ( x ') s «n [t(x) - r(x' )]k £ ( |t(x) - r(x*)|) dr(x') (1.36) 

— OB 

where y is the Stefan-Boltzinann constant and 

, x 

T(x) ‘ J 7 ■ < Wx *I U-37) 

o avg^ ' 

The radiation i3 characterized by the average penetration length. 
This is a consequence of the grey gas assumption. When frequency variation 
of the radiation is considered, there may be no one length which typifies 
radiation energy transfer. Non-grey radiation energy transfer is con- 
sidered in Chapter II. 

Research on radiating shock waves may be classified by whether 
the radiation is lost or trapped. 

Of course no radiation energy can ever be lost in the sense that 
it is destroyed. By lost radiation is meant radiation which escapes the 
system being considered. For example, radiation may be iost from a 
shock tube to the surrounding laboratory. Loss of radiation energy from 
a system is referred to as radiation cooling. 

Radiation which is not lost but re-absorbed within the shock will 
be termed trapped radiation. All previous investigations of shock wave 
structure have considered either trapped radiation or radiation cooling 
but not both. 

Heaslet and Baldwin [6] have obtained solutions for shock waves 
with trapped radiation. The shock wave is assumed to be inviscid and 
propagating in a perfect gas. The mathematical problem associated with 
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the Integral fora of tha radiation flux, which requires the solution 
to be kaovn before the integrand in equation (1.36) can be calculated, 
is avoided by what is essentially the differential approximation 

discussed by Murty [ 7 ] and by Goulard and fraugott [8], 

Although the problem solved by Heaslet and Baldwin is the least 
complicated formulation possible which involves trapped radiation, it 
must still be solved numerically. 

If there were no radiation, the solution would be a step function 
shock wave caused by viscosity and heat conduction in the gas. The nature 
of the solution depends on the amount of heating which occurs by radia- 
tion energy transfer compared to that which occurs by viscous dissipa- 
tion. 

If radiative heating exceeds viscous heating^ it is possible to 
have a shocK wave without any discontinuities. The absence of discon- 
tinuities shows the viscous effects to be completely negligible . Repre- 
sentative profiles of temperature, velocity and radiated heat flux 
typical of this kind of shock wave sketched from the solutions obtained 
by Heaslet and Baldwin are shown in Figure 1,5. Presentation of the 
results in terms of the optical length can be misleading. If the 
average penetration length for grey radiation is large, the extent of 
the shock in terms of the optical length will correspond to a very small 
distance in terms of the physical length x. 

When the radiative heating is not so strong the solution consists 
of a step function shock imbedded in a larger inviscid region. Profiles 
typical of this situation, sketched from the solutions of Heaslet and 
Baldwin, are shown in Figure lA. 
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Heaslet and Baldwin assume the viscous portion of the shock in 
the latter case to be vanishingly thin. Cohen and Clarke [ 9 ] and Chow [10] 
have shown it is a good approximation to consider the shock wave composed 
of an inviscid region, identical to the solution obtained by Heaslet and 
Baldwin, in which is imbedded a viscous region, unaffected b; radia- 
tion, which matches the Inviscid solution at its boundaries. This is 
true as long as the penetration length is large compared to the extent 
of the viscous region. 

Solutions for shock waves with combined radiation and viscosity 
in which there is no restriction placed on the size of the penetration 
length compared to the size of the viscous region have been obtained by 
Traugott [11]. Sen and Guess [12] also solve the combined problem but 
assume the penetration length to be small compared to the thickness of 
the viscous region. They then employ the Rosseland approximation to 
calculate the radiative heat flux. 

All of the solutions just mentioned assume the gas is not ionized 
by the shock wave. Non-grey trapped radiation energy transfer in shock 
waves with ionization has been investigated by Clarke and Ferrari (5l« 
Their work is discussed in Chapter 11. 

Radiation cooling effects in shock tube experiments have been 
reported by Petsch?k et al.[13] and Redkobordyi [14]. 

All of the solutions for radiation cooling assume the shock may 
be divided into two regions. In the first region ionizational equili- 
brium is reached and no radiation cooling occurs. That region is followed 
by the radiation cooling region from which radiation energy is lost. The 
radiation cooling region is terminated by a cold gas interface or wall. 
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Petschek et al. [ IJ] calculate profiles in the radiation cooling 
zone assuming it is optically thin. McChesney and Al-Attar [16] perform 
a similar calculation except they assume in addition the gas is grey. 

in the solution of Pomerantz [ 17 ] the gas is assumed to be grey 
but self absorption is taken into account. In contrast to the optically 
thin solutions, in which the gas continually cools, cooling stops after 
a few penetration lengths due to self absorption. A steady state is 
reached in the radiation cooling region when it is a few penetration 
lengths thick. 

In the case investigated by Pomerantz the radiation cooling region 
was terminated by a cold gas interface. Yoshikawa and Chapman [18] ob- 
tained a solution similar to that of Pomerantz but they terminate the 
radiation cooling region with a wall. In addition to the radiation 
which is lost by emission ahead of the shock there is a loss at the wall. 
Cooling adjacent to the downstream boundary appears to be absent from the 
solution of Pomerantz. 


Chemistry of Ionization 

For shock waves in argon, ionization proceeds in a very complicated 
way. The details of the ionization process are not clearly known. 

Petschek and Byron [19] show that more than one reaction is 
necessary to reach equilibrium ionization. If the electron concentration 
is sufficiently high, the most probable ionization process is ionization 
by electron-atom collisions. Petschek and Byron concentrated on investi- 
gating this reaction. 

They deduced the electron concentration during the approach to 
equilibrium ionization in argon from measurements with potential probes 
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placed In their shock tube and from measurements of the continuum radia- 
tion intensity. Initially the electron concentration is too low for 
electron-atom collisions to be important, but later in the ionization 
process the number of electrons is sufficient for electrom-atom collisions 
to be the main source of ionization. 

Their measured ionization rate is consistent with a two-step 
reaction scheme consisting of excitation followed by ionization of the 
excited state. 

* 

e + A •* A + e 

. . d-5 8 ) 

e+A A +e + e 

Furthermore, their measurements verify that the electrons are not in 
thermal equilibrium with the atoms as ionization proceeds. This is 
because each time an electron is produced the electrons collectively lose 
energy equal to the ionization potential of the atom. The only way in 
which this energy can be regained is through elastic collisions with the 
ions and atoms. However, this is an inefficient process because the 
electron mass is much smaller than the atom mass. On the other hand, 
energy exchange between the ions and atoms, which have essentially the 
same mass, is effective. As a result the electrons tend to cool, but 
the atoms and ions establish an equal and hotter temperature. Petschek 
and Byron assumed the electrons established thermodynamic equilibrium 
qmn n g themselves at an electron temperature T e and the atoms and ions 
established thermodynamic equilibrium among themselves at a temperature 
T . The electron temperature was assumed to be such that the energy 
lost by electrons through ionizing collisions was just balanced by that 
gained from elastic collisions with the atoms and ions. Since the reaction 
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rate la sensitive to the electron temperature, the agreesient between 
measured and calculated rates obtained by Petschek and Byron Justified 
their model of thermal non-equilibrium between the electrons and heavy 
particles. 

Wong [21] has also deduced the electron-atom ionization rate in 
argon by measuring the electron concentration in the chemical non- 
equilibrium region of a shock wave using an interferometric technique. 

He found the electron-atom ionization rate to be more rapid than had 
been measured by Petschek and Byron. He attributes the larger rate to 
the onset of one-step electron-atom ionization 

e + A •* A +e + e 

but it could also be caused by other processes to be discussed presently. 

In the early stages of ionization, when the electron concentra- 
tion is too small to support ionization by electron-atom collisions, 
other reactions must be responsible for ionization. Petschek and Byron 
noticed that the time to reach equilibrium ionization was greatly af- 
fected by the concentration of impurities present in their shock tube. 
Accordingly, they suggested that when the electron concentration was 
small reactions between impurity molecules and argon atoms would be 
important. However, the concentration of impurities in their experi- 
ment was not great enough to account for all of the ionization which 
takes place in the region of small electron concentration. They 
speculated that reactions involving atom-atom collisions and radiation 
in addition to contaminant reactions contributed to the ionization. 

Weymann [21] considered the relative importance of contaminant 


reactions and atom-atom reactions. He shoved the contaminant reactions 
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could be considered Independently from reactions between gas atoms 
because of the much lover concentration of impurities. The impurity 
reactions were conjectured to occur much more rapidly than atom-atom 
reactions due to the lover ionization potential of impurity species. 

He viewed the initial region of ionization as being divided into two 
regions. In the first region ionization would take place by reactions 
involving impurity molecules. This would be followed by a much more 
extensive region where ionization would proceed by reactions involving 
only gas atoms. 

Weymann theorized the most probable atom-atom reaction would be 
a two-step reaction similar to the two-step reaction investigated by 
Petschek and Byron. 

A + A 
A* + A 

The concentration of excited itoms is maintained either because the 
excited state is metastable or because radiative transitions from this 
excited state are optically thick. This reaction scheme was experimental- 
ly verified by Harwell and Jahn [22] by use of a microwave probe trans- 
verse to their shock tube. Great care was taken to reduce the impurity 
concentration. Only when the impurity level was reduced to a few parts 
per million did their presence have a negligible effect on the results. 
Their technique for reducing the impurity level was further refined by 
Kelly [25] who confirmed the conclusions of Harwell and Jahn and obtained 


* 

A + A 


<* A + e + A 


(1.39) 


more accurate results. 
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The contaminant reaction region was lnveatlgated In xenon by 
Backer and Bloomberg [2k] by using the microwave probe along the axis 
of their shock tube. They found the ionization occurred by a very 
complicated set of reactions involving contaminant species, xenon atoms 
and molecules in various states of excitation, and radiation. Contami- 
nant reactions in argon were theoretically discussed by Morgan and 
Morrison [25] who calculated the time to reach equilibrium ionization in 
argon shock waves by combining the electron-atom rate of Petschek and 
Byron and the atom-atom rate of Harwell and Jahn. They attributed the 
discrepancy between the calculated time and the time observed by Petschek 
and Byron to contaminant reactions. The calculations made in this thesis 
assume the argon gas to be absolutely free of impurities. 

Blbermann and Yakubov [26] also calculated the time to reach 
equilibrium ionization by combining electron-atom ionization and atom- 
atom ionization rates. Experiments show the ionization process actually 
takes place more quickly than they calculated. They suggested that line 
radiation from the hot gas following the ionization region excites argon 
atoms in the ionization region, which are easily ionized, contributing 
additional electrons and consequently decreasing the time to reach 
equilibrium. The equilibrium ionization time calculated by including 
this effect showed much improved agreement with experimentally measured 
times. 

Radiation can affect ionization in still another way by creating 
electrons ahead of the shock discontinuity. Such electrons have been 
termed precursor electrons because they occur ahead of what is normally 
thought of as the shock front. They have been observed in re-entry by 
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their effect on redar croee section [27], (28], [29], and In shock 
tubes as described belov. 

Results of experiments to detect precursor electrons In shock 
tube6 were at first quite confusing. Some care was needed to eliminate 
extraneous causes of electrons ahead of the shock wave such as the photo- 
emission from the shock tube walls [30], [3-1]* The early measurements 
of precursor phenomena seemed to indicate that most of the charged parti- 
cles ahead of the shock front were electrons [30], [32], [33]. If 
the photoionisation caused the precursors^ the number of electrons and 
ions would be equal. It was proposed that precursors were caused by 
electrons diffusing ahead of the shock wave. It was thought the severe 
gradients of electron concentration in the shock wave would cause such a 
diffusion. 

Theoretical investigation of the diffusion hypothesis has been 
reviewed by Wetzel [3^]. AH the calculations snowed a much soiailer 
electron density ahead of the shock than had been measured. Pipkin {35] 
argued that the distribution of electrons in small diameter shock tubes, 
such as had been used for the experiments, produce an electric field 
vector that Is not mainly along the axis of the shock tube as had been 
assumed in the theoretical work. He analyzed the problem assuming a 
different configuration of the electric field and, because of its weaker 
axial component, was able to explain the main features of the experi- 
mental results. Appleton [36] claims it is unreasonable to assume that 
Pipkin's expression for the axial field is valid in the immediate vicinity 
of the shock. Instead Appleton requires the field to be one-dimensional 
near the shock bnd obtains a numerical solution in which the electron 
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concentration rapidly decreaaea ahead of the abode, conaiatent with the 
earlier theoretical result*. Unless a transition to a weak axial field 
takes place in a short distance ahead of the shock, the diffusion hypothe- 
sis cannot be theoretically explained. 

Recently Weyoann ana Holmes [37] and Holmes [ 36 ] have repeated 
the earlier experiments of Weyaann with Improved measurements. These 
more recent experiments show that photoionization of the cold gas ahead 
of the shock by radiation from the hot shocked gas is the principal means 
of producing precursor electrons. This is consistent with the recent 
experiment of Zivmaovlc [ 39 ] an a Lederman and Wilson [kOj who were able 
to separately wea^m e the photo-electrons and the diffusion electrons 
ahead of the shock, waves. They found only photo-electrons. 

These experiments indicate that precursor ionization ia caused by 
radiation from the hot shocked gas. 
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CHAPTER II 

FORMUIATION OF EQUATIONS 

In this chapter the equations which describe non-grey radiation 
energy transfer, ionization kinetics, and energy exchange in partially 
Ionized argon are formulated. First the fluid conservation equations 
for a partially ionized gas are developed. These are then transformed to 
a non-dimensional form and simplified by eliminating all unimportant 
terms. Next equations for the ionization rate are developed. Finally 
expressions are developed which give the rx>n-grey radiation flux and 
radiative ionization rate. The method of solution is discussed at the 
end of this chapter. 

Conservation Equations for a Partially Ionized Gas 
In an un-ionlzed monatomic gas mass, momentum and energy can be 
transferred only by the neutral atoms which make up the gas. If the gas 
is ionized, electrons and ions provide additional ways for mass, momentum 
and energy to be exchanged. An equation for the conservation of mass, 
momentum and energy can be written for each species of the gas following 
the approach used by Appleton and Bray [4lJ. Let each species be labeled 
by a subscript X. The species conservation equations for one -dimensional 
steady state flow are as follows: 

The species mass conservation equation is 

If < »n V ■ *x 


(£. 1 ) 
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where ■ m^n^ it the aass density of epeciet X, m^ Is the mass of 
a particle of species X, is the number density of species X, is 
the average velocity of species X, and is the mass rate of production 
of species X per unit voluee per unit time. 

\ 

where is the net number of species X produced per unit voluae per 
unit time by all radiative and collisions! reactions. 

The equation expressing conservation of momentum for species X is 



n. 


e x E 


l 

x» 


+ J, 


XX' 


( 2 . 2 ) 


where p^ is the partial pressure of species X, is the viscous stress 
tensor component of species X, e^ is the charge of a particle of species X, 
B is the electric field created by the collective motion of the charged 
particles, P is the rate of lost, of momentum by species X per unit 

I 

volume and time due to elastic interactions with particles of species X, 
and is the rate of loss of momentum by species X per unit volume and 
time due to inelastic interactions with particles of all other species 
and the radiation field. 

Conservation of energy for species \ is expressed by the following 
equation. 


s'l'r i p x u x * - °x> u x*\> 

• n x e x Eu x * u x ^ P XX. + I V + S 


where e is the average internal energy per unit volume due to species X 
A 

and q c is the conductive beat flux. The first term on the right hand 
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side of equation (2.3) is the energy dissipated by Joule heating and the 
second is the rate at which work is done on species X per unit volune and 
time due to elastic interactions between species X and the other species be- 
cause of the fluid motion of species X. In the remaining terms E^, is the 
rate of energy gain by species X per unit volume and time due to elastic 
encounters between species X and X* because of thermal motion of the parti- 
cles and is the rate of energy gain by species X per unit volume and time 
due to inelastic encounters between species X and all other species because 
of the thermal motion of the particles and radiation. 

Equations (2.1) and (2.2) may be used to reduce equation (2.3) to 
the following alternative form 


d_ 

dx 


(e xV \ ] + 





**XX' 



U.*) 


Elastic Interactions 

Assuming that the degree of ionisation in the gas is sufficiently 
small that only binary collisions are important, the following kinetic 
theo y expressions give the net momentum and energy transfer between species 


X and x* [lu], [42]. 


P XX' = " ny W j*Kx ,{g) f X f V dv X dV X' 


^XX' " “ Vx'J S XX ,(s) f X f X’ dV X dV X’ " Vxx* 


(2.5) 


( 2 . 6 ) 


In these expressions g is the relative velocity of the particles before 
collision, G is the center of mass velocity of the particles, is tne 
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velocity distribution function of species X, is the total velocity 
of species X equal to u^ + c^ where is the thermal velocity, and 
S^, is the diffusion cross section between species X and X’. 

It is reasonable to assume the time required for each species 
to establish a Maxwellian distribution of thermal velocities by self 
collisions is short compared to the time required to exchange energy with 
other particles [45]. Accordingly, the velocity distribution function 
in equations (2.5) and (2.6) are assumed to be Maxwellian at a particular 


species temperature T, . The species temperatures are equal only if 

thermodynamic equilibrium between the species is established. 

In addition to the distribution functions, the cross sections 

S^.., are needed to evaluate the integrals in equations (2.5) and (2.6). 
AA f 

They are given by. the following equations [4], [42], [44] -[47]. 
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where the subscripts 1, a, and e stand for ions, atoms, and electrons. 

Using these cross sections and the Maxwellian velocity distribu- 
tion functions the integrals in equations (2.5) and (2.6), correct to 
first order in (m /m )^, are [fc], [ 19 ]# l-^8],l5l3» 
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Pressure, Viscosity and Heat Conduction 
The partial pressure viscous Btress and conductive flux q Q 
are assumed to be given by the following familiar formulas. 


»x * “x k \ 


( 2 . 18 ) 
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X * 3 Wj x dx 


(2.19) 


dT. 

q - - k ~ 

C X C X 


( 2 . 20 ) 


The viscosity coefficient ii and thermal conductivity coefficient 

U X 

k are given by the following classical expressions [4], [48). 

C X 
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( 2 . 22 ) 


In these equations c is the mean thermal velocity of species X, 


\ J ^ ) 


(2.23) 


and is the frequency of collisions of species V with all other 

particles. 
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Noting that m /m « 1, the viscosity and heat conduction coef- 
6 & 


ficients may be written as follows: 
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The relative magnitude of these coefficients depends on the amount 
of ioniaation that has taken place through the ratios n.s. /n P. 

^ 1 lA 9 98 ? 

n i s ii /n a 3 ia> etc ‘ 


Inelastic Interactions 

Inelastic interactions are collisions in which the internal energy 
of the collision partners is permanently changed, s’ach species may 
gain or lose energy through inelastic collision processes. 

Processes of this kind are the inelastic electron-atom (e-A) colli- 
sion process 

e + A(l) Z e + A(2) 
e + A(2) Z A + + e + e 

and the inelastic atom-atom (A-A) collision process. 

A (1) + A(l) Z A (2 ) + A( 1) 

A(l) + A(2) Z A + + e + A( It 

Particles may also exchange energy with the radiation field. 

Consider how the e-A process will affect the particles in a unit 

volume of gas. There are n electrons in the unit volume. Suppose one 
0 e 

electron is produced by the e -A process. The electrons which cause the 


52 


Ionization will lose energy equal to the ionization potential X plus 

the energy of the created electron. The average energy of the created 

electron is The original n electrons will establish equilibrium 

by elastic collisions among themselves. Thus the n original electrons 

© 

collectively lose energy X + £ in creating one electron of average 
energy There will then be n 4 1 electrons in the unit volume, 

n Q of which have lost the energy )t + f and one which has gained energy 

£. No energy is lost in the unit volume by equilibrating the created 

* _ 

electron since the required energy ^ k - £ will be lost by the 
original n^ electrons but gained by the created electron. The net 
energy loss per unit volume is therefore Just X. If n g ^ is the net 
rate of production of electrons by the e-A process, the total energy loss 
per unit volume by this means is * n eA *» 

Since 1 k T 4 u' is the average energy of an atom, and 

the e-A process destroys atoms at the rate -n^, the atoms will collec- 
tively lose energy at the rate -n ( i k T 4 x m uf), Similarly 

68 C 8 c. o. Ql 

the ion gas will gain energy at the rate n ga ( ^ k ^ m^u^). 

Electrons are Ineffective in transferring momentum with ions and 

atoms because of their small mass. Approximately the rate of momentum 

exchange by the e-A process for the electrons is zero, for the atoms is 

-n . m u and for the ions is n . m. u . . 
eA a a eA i i 

Similar considerations for the A-A process and radiation lead to 
the following formulas for the inelastic momentum and energy exchange. 
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where £ is the average energy of electrons when created by the A-A 

process, is the average creation energy of photoelectrons, n^ is the 

rate of production of electrons by the A-A process, n ^ is the rate of 

production of electrons by the radiation field and n«n.+n AA +n , . 

eA AA rad 


The Electric Field 

Steep gradients of electron and ion densities in the shock wave 
will cause diffusion of electrons relative to the ions. The resulting 
charge separation creates an electric field. The electric field is found 
from Poisson’s equation; 


dS 

dx 


( n i ' il e ) 


( 2 . 57 ) 


The diffusion is impeded by the electric field it creates. A steady 
state is reached when the diffusion is Just balanced by the drift from the 
electric field. It is assumed a steady state diffusion is maintained in the 
shock wave. The net current then is zero. 
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Internal Energy 

The Internal energy due to the thermal motion of particles of species 
X Is 3/2kT^ since the particles have a Maxwellian velocity distribution 


at the temperature T^. The internal energies of the atoms, 
trons per unit volume are: 

ions and elec 


’ “a 

(2.59) 
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where ^ is the average electronic excitation energy of the atoms. For 
the calculations to be performed t|r ft is negligible compared to the 
thermal energy of the atom. 

Simplified Nondimensional Equations 
It will be advantageous to refer the variables to characteristic 
values in order to aid in the estimation of the relative magnitudes of 
terms in the equations. Let the characteristic values be denoted by a 
subscript U. Then the barred variables 



will all have a value near one. Let l be a length which is typical of 
the extent of the chemical and radiative portion of the shock wave. It 
will be on the order of a centimeter for the shock waves to be considered 
in this thesis. Let x = x ,- 1 be the nondimensional coordinate. 

In addition to the length l there are several other characteristic 
lengths which will appear in the equations. Those which are associated 
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with the momentum and energy exchange are presented next, 

in the course of the presentation to follow it will be shown that 
charge separation is negligible for the shock waves to be investigated. 
Then a • n. and the species concentration may be related to the decree 
of ionization by the following formulas. 
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Using relations (2.4?) the lengths characteristic of the energy and 


momentum exchange processes may be written as follows. 
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(l - a)p a „(£*). 


(2.50) 
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and S ea (T e ) ■ S^ a (T e )/fl^ a (T 0 ). Using the rsference values T Q - 10* *K 
and a density of p Q ■ 2.15 (10 ), gm/ cm (which is the density corre- 

sponding to a pressure of 1 cm Hg at a temperature of 300 *K) the following 


values are found. 


1° » 1.31 (10‘ 4 ) cm 

SB 

= 1.58 (IQ'S cm 


14.5 cm 
.8 (10* 2 ) cm 


(2.52) 


From these formulas it can be seen that the lengths for heavy 
particle partners, i and l^ a , are considerably smaller than those 
involving electrons, / #a and . 

An appropriate nondimensional electric field is 


(4 V e 


(2.53) 


where ijj is the Debye length. In terms of nondimensional variables 


the Debye length is 


where 
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For the same values of the reference variable as used above, 1° la 
found to be 1.21 (10*^) cm. The Debye length la auch smaller than any 
other length in the problem. 

Sutoattttttlig the f omuls* developed In the last section and the 
nondlmensional variables from this section Into the electron momentum 
equation (equation (2.2) with X » e) and noting that m « m produces 
the following equation i 


d__ 

dx 
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where the parameter a « kTQ/m a UQ is a typical ratio of thermal to 
kinetic energy of the particles. Since Ip « I, the first term on 
the right hand side will be the dominant term in the equation. Unless 
the gradients in n a ml T are large the electric field must be small in 
order to preserve the equality. It will be assumed that S is small 
throughout the region which has characteristic length I. Equations 
(2.37) and (2.3 d) then show that there is no charge separation and no 
diffusion of electrons relative to ions. 


\ - n e (2.57) 

\ ( 2 . 58 ) 

This justifies the previous use of relations (2, 4i) which will be used 
frequently in the equations which follow. 

Consider now the equation obtained by adding the electron momentum 
equation to the ion momentum equation. This is found by summing equation 
(2.2) over the electron and ion species. Using the result just obtained, 
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n. ■ n and u. • u, , and tha formulas specified in the previous sections, 
the nondi— naional form of this equation is found to be 


a "7 (5 %) ♦ ~ I ' • P 

ax dx L 



Since the ratio -r — , which appears on the right hand side of the 

xa 

equation, is very large, the mean atom velocity u a must be very nearly 
the same as the mean electron velocity in order that the equality 

be true. Consequently there is negligible diffusion of the species 
relative to each other. 
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The largest term on the left hand side of this equation is of 

order a. On the right hand, side of the equation the term containing 

£ . 

the ratio 7- and the term containing the chemical reaction rate n 

^ia 

will have the largest values. If the chemical rate were zero, the atom 

temperature and the ion temperature would have to be very nearly equal 

to preserve the equality since the ratio 7 — is very large. However. 

*ia 

if the ch emical rate is sufficiently high, the atom and ion temperature, 
can iffer. For the shock waves which are considered in this thesis the 


chemical rates are sufficiently small that differences in the atom and 
ion temperatures are negligible. Accordingly, it is assumed that 


T a = T ± (2.62) 

The assumptions just discussed greatly simplify the problem for 
they have reducea the dependent variables to the following five quanti- 
ties: a, p, u, T , and Te. The equations which will be used to calcu- 

late these variables will be the total conservation equations, the 
electron energy equation and electron mass conservation equation. 
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The total conservation equations are obtained by sunning equa- 
tions (2.1), (2.2) and (2.?) over all species. 

The noadimensional fora of the total mass conservation equation 
is 

P U ■ Cjj m 1 (2.65) 


The constant Is chosen to be one. This means that the reference 
quantities must be chosen such that p q u q is equal to pu at some point 
in the shock wave. 

After summing equation (2.2) over all species, substituting the 
nondimensional variables and integrating once, the following equation 
which governs the transfer of momentum is found. 
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where C p is the integration constant. The last term in this equation 

*aa 

is the viscous stress. This term is negligible since the ratios r— 

I . * 

and rr-~ are very small and the velocity gradient is assumed not 

* dx 

to be large. The equation of momentum transfer for the gas is therefore 


—— + a p (T a + a T e ) * C p (2.6b) 

The equation of total energy conservation is obtained in a manner 
similar to that used to derive the total momentum conservation equation. 


The result is 
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(2.65) 


where Cg is a constant resulting from the integration of equation 
(2.5). The last two terms represent thermal heat conduction in the 
gas and the one preceding these represents viscous dissipation. The 

* aa *i a 

terms containing the ratios — — and p— are clearly negligible but 
the term containing the gradient in electron temperature car be impor- 
tant. This term represents heat conduction in the electron gas. Its 
larger magnitude reflects -he fact that the electrons are good con- 
ductors of heat. As can be seen from equation (2.47), the quantity 
J, 


a 


ee 


, which is the magnitude of this term, is sufficiently small for 


the range of variables of interest in this thesis that heat conduction 
in '.he electron gas may be neglected. The equation of energy conserva- 
tion then reduces to the following form. 
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(2.66) 
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Equation (2.4) evaluated for X » e is the equation which de- 
scribes the energy balance for the electrons. Substituting tb*. 
formulas developed in the preceding sections into this equation and 
noting that n./a « 1. the following nondiaensional fora of the elec* 
tron energy equation will be found. 



o o 


Conduction and convection of energy in the electron gas is 
described by the terms on the left hand side of this equation. The 
gain in energy of the electrons by elastic interactions, inelastic 
interactions, and the creation of electrons by ionizing reactions is 
accounted for by the terms on the right band side of the equation. 

When the rate of production of electrons is sufficiently slow 
the right hand side of equation (2.67) reduces to the single term 

This terra represents the energy gained by electrons because of elastic 
interactions between the electrons, ions and atoms. 


L 
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Jaffrin [ 4 ] has studied shock waves in ionised argon assuming the 

rate of production of electrons to be zero. Equation ( 2 . 67 ) reduces to 

the electron energy equation used by Jaffrin when the rate of production 

of electrons is zero. In addition Jaffrin neglects the radiative heat 

flux q^. He finds the shock wave consists of a region of the order 

i thick in which the electrons are heated by conduction and elastic 

interactions with ions and atoms. Imbedded in the region where heating 

of the electrons takes place is a much smaller region of the order l la 

or I thick in which the atoms and ions experience a viscous shock. 

For the shock waves studied in this thesis the rate of 

production of electrons is not negligibly small. Heating of the 

electrons by conduction and convection, however, is negligible because 

^ee 

either a is small, or a •j— is small or the gradients in T g and u 
are small. Consequently the electron energy equation reduces to 
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( 2 . 68 ) 


When all ionization processes except the electron atom ioniza- 
tion process are ignored equation (2.68) reduces to the equation used 
by Petschek and Byron [19] to describe the energy balance of the 
electrons . 


The final equation needed to completely determine the problem 
is the mass conservation equation for the electrons (equation (2.1) 
with "K ■ e). The nondimensional form of this equation is 
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Sunnary of Equations 


The simplified nondimensional equations which are used in this 
thesis to describe shock waves in radiating and argon are: 

The total mass conservation equation. 

p u * 1 (2.6: 

The total momentum conservation equation. 


+ a p (T # + a T e ) - C p 


The total energy conservation equation. 
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The electron energy conservation equation. 
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The electron mass conservation equation. 


da n 


It remains to give detailed expressions for the ionisation rates 
and the radiative heat flux. This is done in the following two sections 


Ionisation Rates 


>»5 


In order to compute the rate 8 h^, n^ and n^^ it is 
necessary to know the population of all electronic energy Avals' of the 
atom which contribute to the rates. In many instances one process, for 
example, electron collisions, is dominant in populating these states and 
a Boltzmann distribution of the population of excited levels can be 
assumed. 

However, in the present problem no one process is necessarily 
dominant in populating all the electronic states of the atom. For 
example the upper levels might be affected most by collisions of the atom 
with electrons, but the ground level population might be affected most 
by radiative interactions. 

In theory it is necessary to calculate the population distribu- 
tion of the atom by solving a rate equation for each level, taking into 
account all the ways by which the population of that level can be 
changed. Such a calculation has been performed for atomic hydrogen 
(the only gas for which the collisional and radiative cross sections are 
known in sufficient detail) assuming the upper states of the atom to be 
in statistical equilibrium [52} - [5 6], 

A detailed analysis of the population of excited levels in 
argon cannot be made because of an insufficient knowledge of the neces- 
sary data. Kven if all the necessary data were available it would be 
desirable to simplify the calculation by making a reasonable model of 
the actual case. 

The model employed in this thesis is based on the fact that 
electron-atom collision rates greatly exceed other rates for those 


levels which are close to the ionisation Unit of the atom [52]- [ 58 ]* 
Consequently these levels will establish equilibrium with the electrons. 
Since the excited states of the argon atom are all within about 4 ev of 
the ionisation limit, the excited states of the argon atom are assumed 
to be in thermodynamic equilibrium with the electrons. 

Since the excited states are in equilibrium with the electrons, 
their net rate of production of electrons is zero. All of the free 
electrons must be supplied by the ground level. Ground state electrons 
are released by electron-atom collisions, atom-atom collisions and 
radiative interactions. 


The colli slonal reactions 

depleting the ground state are 



K A (l,g) 



A(l) + A 


A(2) + A 

(2.70) 





A(l) + e 

K.(2,D 

A(2) + e 

(2.71) 


where index 1 denotes the ground state and index 2 the excited state. 
These reactions are followed by a rapid ionization from the excited 
state. 

The rate of change of population of state 2 due to atom-atom 
collisions is 

A (2) = K a ( 1,2) n a n (1) - ^(2,1) n(2^ n a - n + (2) 

where n + (2) is the rate of ionization from state 2. Since state 2 is 
in equilibrium with the electrons, its net ionization rate must be zero. 
Consequently the rate of ionization by atom-atom collisions is 


A AA • V 1 ' 2 ) \ n < 2 ) \ 

where the rate coefficient K^(l,2) is shown in Appendix A to be 

2C* 3/2 / T r ._ \ T ex/ T a 


(2.72) 
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(2.73) 


In equation (2.73) T ex ■ 1, 34(l(r ) *K is the excitation temperature 
for the first excited state of the argon atom and ■ 7.5(10 )cm /erg 
is a constant used in the expression for the atom-atom excitation 
cross section. 

At equilioriura must be zero and the population of the 
atomic levels must be given by a Boltzmann distribution at the atom 
temperature T . In tne limit of equilibrium equation (2.73) ehows that 

A 

K^(l,2) is related to K^(2.1) by 
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where is the degeneracy of state 1. Consequently the atom-atom 
ionization rate may be written as 


*AA ‘ M 1 . 2 ) V (1 > 1 1 1 

“ ‘ L ln(2)/n(l)J u « J 
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(2.75) 


Since state 2 is in equilibrium with the electrons, its population 


may be found from tne Saha equation. 
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where it the ionization potential of state 2. 

The population of state 1 which would exist if it were in 

equilibrium with the electrons n(l) «. , is found from the follow- 

b, 

ing similar equation 

2 5/2 

2(2 n m e k T.) - X/kT e 


n(I) 


S, T e 


g A h- 


Therefore the ratio of the population of states 1 and 2 may be written as 

SlSi * n( ^SA h 'V 1 . 

SfiT g x e 

The firs,; excited level of the argon atom is separated from the ground 
state by an energy gap of about 11.8 ev. Because of this large energy 
gap, it is a good approximation tiiat 


n(l) - « (1 - a) o/m ft 


(2.76) 
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(2.77) 


where Ctg(T r J is the degree of ionization which would exist if the gas 
were in equilibrium at the electron temperature and is given by the 
following formula. 
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(2.78) 




The ratio of the population of states 1 and 2 aay then be written as 


m . .- T « /T « s (2.79) 

1 ' l-o *1 

Using equation (9.79) in equation (2.79) the atom-atom ionisation rate ie 
V „ N2 r 1 - o (T.) - T ex(-?; ’ -A 

°*A ‘ (1 ' ®) L 1 ‘ 1 ~ - ' J ~ * J (2.80) 


In a similar way the electron-atom ionisation rate is 
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( 2 . 81 ) 


where K e (l,2) is the forward reaction rate constant for the electron- 
atom reaction. It is shown in Appendix A to be 
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( 2 . 82 ) 
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In equation (2.82) C g * 4.4 (10” ) cm^/erg is a constant used in the 
formula for the electron-atom excitation cross section. 

When the degree of ionization is small the atom-atom rate is 
larger than the electron-atom rate. By substituting equation (2.60) 
into equation (2.66) and assuming the degree of ionization is small, 
it can be seen that as the degree of ionization decreases the electron 
temperature also decreases. Consequently for small degrees of ioniza- 
tion the collisional ionization rate reduces to 


<V • U - a) 2 ( *■ ) *(1.2) 


It may be that just behind the viscous shock the degree of 


ionization is so small that there are insufficient electrons to establish 
equilibrium with the excited states of the atom. Equation (2.79) 
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then invalid. Just behind the viscous shock, however, it will be true 
that 

« 1 

(n(2)/n(2)) 

» a 

and equation (2.75) shows that the colli sional rate will be the sane as 
the expression just given for low degrees of ionization, therefore 
although equations ( 2 . 60 ) and (2.61) are invalid for very small degrees 
of ionization, their limiting values are correct. Consequently the ioni- 
zation rates will be assumed to be those given by equations ( 2 . 60 ) and 
(2.61) for all degrees of ionization. 

In the equations developed in the last section the ionisation 
rate always appears divided by the characteristic convection rate N. 
Dividing equations (2.60) and (2.61) by N and expressing the result in 
terms of nondimensional variables gives the following equations. 
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■otic* that the electron-atom rate depends on the electron 
temperature and the atom-atom rate depends primarily on the atom tempera- 
ture. When equilibrium is reached T e » T ft and a • a^(T e ). The rates 
then are zero. 

The radiation ionization rate is derived after the next section. 

Mon-Grey Radiation Energy Transfer 
The ways in which radiative interactions take place in a non- 

S 

atomic gas are described by the following reaction equations. 

A(p) ♦ h v A* + e p = 1,2,3 (2.86) 

and 

A(p) ♦ h v A(q) q > P p = 1,2,3 (2.87) 

where A(p) denotes an atom excited to the electronic state. 

Appendix B shows it is possible to write a radiation transfer 
equation for each of the radiative reactions in equations (2.86) and 
(2.87) and each radiation transfer equation has the same form as equation 
(l.2k). Let the radiative reactions be labeled by the subscript j. The 
radiation transfer equation which corresponds to reaction j is 

dl. (x,u,v) . 

■•H 4 ’ ' 7^x"v7 [I(X ' * v) - s j (x - v)1 (2 - 88) 

The manning of the terms in this equation is exactly the same as for 
the corresponding terms in equation (1.24) except they apply to the 
j*’* 1 reaction only. For example I (x, u, v) is the intensity of radla- 

J 

tion from the Jth reaction alone. 
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The total intensity Is the bub of the individual intensities 
fro* ell the radiative reactions. 


I(x # n, v) - 


(2.89) 


If equation (2.86) is sunned over ell the radiative reactions, equation 


(1.24) is obtained. 

Z AIj(x, n, v) 


»** v) 


* u 


i (x,v) 


U(x,H, V) - S(x, v)] 


where the penetration length l(x,v) and the source function S(x,v) are 

the following functions of the penetration lengths l^(x,v) and source 

functions ?>.(x,v) for the individual radiative reactions. 

0 


1 _ V j- 

(x,vT " L "Tlx/ 


s(x ' v) • s j ( x ' v) 


(2.90) 


(2.91) 


Once the individual penetration lengths l ,(x,v) and source functions 
S (x,v) have been specified, equations (2.9°) and (2.9I) can be used in 
equations (1.26) and (1.52) to calculate the radiative beat flux. 

For the photoionization reactions (2.86) the following formulas 
for the individual penetration lengths and source functions are derived 
in Appendix B. 
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(2.95) 


where n(p)g is the population of atoms in state p which would exist 
if level p were in equilibrium with the electrons and ions. It is 
found from the following formula - 

5/2 

- X_/k T 

( 2 . 94 ) 
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p 

In equation (2.<A) ^ is the ionization energy at level p and g^ is 

the degeneracy of state p. 

The quantities B (v) in equation (2.92) represent cross sections 
pc 

for photoionization from level p. In the derivation of equations (2.92) 
and (2.95) the electron energy was assumed to have a Boltzmann distribution 
at temperature T g . 

As in the sections on ionization rates all but the ground state 
are assumed to be in equilibrium with the electrons. 


a(p) B = n(p) 


P ® 2, 5»4, ... 


and for the ground state 
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Using these simplifications the penetration length and source 
function for photoionization from the ground state are found to be 
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The source function given by equation ( 2 t $ 6 ) differs sppreclebly from the 
Planck function only when a g (T ft ) is near one and o / Og(T e ). 

For photoionisation from excited states the following formulae 
are found* 

1 
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* B pc (v) n(p) (1 - e ) P * 2,?,k ••• 


( 2 . 97 ) 


ov 2 hv/kT , 

8 (x,v) - (e 6 - l)* 1 » B v (T e ) p- 2,3,4... (2.98) 

C 

It is to be expected that the source function for photoionization 
from excited states, equation (2.98), is the Planck function since the 
excited state atoms have been assumed to be in equilibrium with the 
electrons. 

The penetration lengths and source functions for the photo ioniza- 
tion reactions ( 2 . 87 ) are shown in Appendix B to be 
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( 2 . 99 ) 


( 2 . 100 ) 


where B is the total transition probability for the transition between 
states p and q, and P (v) is the line shape function. 

pq 


Separation of Radiative Interactions 
Radiation energy transfer is combined in the sense that the total 
penetration length and source function depend on the penetration lengths 
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tad source functions for etch of tbs individual radiative interactions 
as shown by equations (2.90) and (2.91). 

However, if for sons range of frequencies one penetration length 
is significantly snaller than all the others, the tern in the stagflation 
containing that penetration length will be dominant. For example if 
J k (x,v) is considerably smaller than all other penetration lengths for 
the frequency range v ft * v S v fa , equations (2. 90) and (2.91) will 
simplify to 

S(x,v) - S k (x,v) 

1 1 
r(x,v) * . j 

These equations show radiation vould be transferred in the frequency 
range v a S v i as if the k^* 1 reaction were the only radiative 
interaction taking place. 

If the penetration length /^(x, v) is assumed to be vanishingly 
small compared to the scale of variation of gas properties, equations 
(1.26) and (1.32) show the radiative flux for the k reaction is zero. 
Consequently reactions for which the penetration lengths are vanishingly 
small may be eliminated from the summations in equations (2*90) and 
( 2 . 91 ). 

Radiative reactions for which the penetration lengths are very 
large may also be neglected since terms for which the penetration lengths 
are large will not contribute to the summations in equations (2.90) and 

( 2 . 91 ). 


V a s V s v b (2.101) 

V a s V S V b (2.102) 
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Line Radiation . The photoexcitation reactions (2*87) product 
line radiation. Many of the Uses may be eliminated from conaideration 
because the atomic levels shram which the excitation proceeds are mot 
populated* Equation (2.99) above the penetration length is then very 
large. 

For those lines produced by photoexcitation from populated states 
the penetration length at the line center is usually vanishingly small 
and radiation energy transfer in the center of the line may be ignored. 
However, the penetration length in the line wings may be such that 
radiation energy transfer is important. 

Murty (15] has examined line radiation in the precursor region 
ahead of shoe* waves in hydrogen. He finds the excitation of atoms 
caused by the line radiation to be important to the ionization process 
far ahead of the shock front. Biberman and Yakubov [26] have concluded 
that radiative excitation of atoms affects the approach to equilibrium 
ionization for shock waves la argon. 

However, it is advantageous to determine the effect of continuum 
radiation without the considerable complication which accompanies the 
analysis of line radiation, For this reason further consideration of 
line radiation is omitted. 

Continuum Radiation . The cross section for photolonisatlon from 
the ground state, B^(v), is zero for frequencies less than the fre- 
quency of the ionization edge, « X/h. Consequently for frequencies 
in the range 0 s v a only the radiation due to photolonisatlon from 
excited atomic states takes place. Equations (2.90), (2.91), (2.97), and 
(2.93) then give 
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S(x,v) ■ B v (T t ) 


0 < v < v x 


(2.103) 


. -hv/kT r-> 

* (1 - • > l V < v) °<*> 

P * 2 

Shock tube experiment e in argon show the continuum radiation ia 
well represented by the Unsold formula [ 13], [14], [371. The basic 
assumption used to obtain the Uns&ld formula is that all the excited 
levels are hydrogenic and closely spaced. Mies [ 38 ] and Dronov et al. 
[59} improve tills model by assuming all the levels to be hydrogenic 
except for certain lower lying states which are treated in more detail. 

In view of the good agreement of the Unsbld formula the refinements of 
Mies and of Dronov et al« need not be used for the shock waves of 
interest in thi6 thesis. 

The Unsbld formula predicts the emitted radiation is constant at 
all frequencies. Since the radiation must decay at the higher frequen- 
cies, the radiation predicted by the Unsbld formula is usually terminated 
at some rather arbitrary cutoff frequency. However, in the formulas 
developed below, which are based on the same assumptions as the Unsold 
formula, the radiation properly decays at the higher frequencies. 

Assuming the excited states are hydrogenic, the Menzel and 
Pekeris formula [60j for the photoionization cross section may be used 
in equation (2.104) to give 
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(2.105) 
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where K ■ ■ ■- y - ■— R^, R^ is the Rydberg constant, is the 

*• ® 5 v 5 

effective nuclear charge, 0 the Gaunt factor, i) the energy of the photo- 
electron. the siamation over the atonic state in equation ( 2 . 105 ) is to 
be carried out at constant frequency • Consequently, the electron energy 
is limited to values consistent with the requirement that the transition 
energy, hv, be a constant 


hv 


K * n 


constant 


(2.106) 


where X is the energy of level p below the ionization limit 
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(2.107) 


Assume the states are closely spaced down to a level X^ below 
the ionization limit and regard the quantum number p as continuously 
distributed in this range. Then from equations ( 2 . 107 ) it follows that 
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eff 
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and since hv * constant * X + q 
dX m . dq 

The maximum electron energy consistent with equation (2.106) is 
that corresponding to X « o and is 


n = hv 

'max 


The minimum allowed electron energy consistent with equation (2.106) la 
that corresponding to the largest value of X which is X^. 
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Substitution of these results into equation (2.103) gives 
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As tuning the Gaunt factor to have a constant value and performing the 
integration the result is 
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■ 1 if V < V M ■ *y|/h 

f { B (T ) v M 3 
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At frequencies lets than the frequency v M the penetration length 
given by formula (2.108) is the tame as that which would be found from 
the Unsold formula. At frequencies greater than the penetration 
length given by equation ( 2 . 108 ) differs from the Unsold result by the 
factor f. Similarly for frequencies less than v^/the intensity of 
emitted radiation, which is B y (T e )/l(x,v), is independent of frequency, 
but for frequencies greater than v M the intensity of emitted radiation 
dscays as l/v 3 . The penetration length has a marked frequency dependence 
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through the Planck function lad tha factor f . 

How consider ground state photoionizatlom* Hie cross section 
for photo ionisation from the ground state of argon has been experi- 
mentally dete mined [61]. It can be roughly represented by 


B l C ( y ) 


V v l J 


if v < » x/h 


if y > « X/h 


where = 35(10” / cm' is the cross section at the ionization 

threshold . 

For frequencies equal to and larger than the penetration 

length for photoionisation from the excited states, equation (£.108), 
will be very large. Consequently excited state pbotoionization 
may be neglected compared to the ground state photoionization 
when v £ v^. 

Collecting the results of this section, the penetration length 
and source function to be used in the calculation of the radiative 
flux are 
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(2.U1) 


Clarke aad Ferrari [ 5] have also studied non-grey radiation 
energy transfer in shockwaves with ionization. Their work differs 
from this study in that they use different reactions to describe the 
collisional ionization process; they assume the electron and atom 
temperatures are equal and they neglect photoionization from excited 
states. 


Reduction to Grey Radiation 

Non-grey radiation energy transfer can be reduced to an equiva- 
lent grey radiation energy transfer if the penetration length can be 
written as a function of x times a function of v. Assuming this 
write 

i(x, v) * 0(x) (2.112) 

The contribution to the radiative heat flux in the frequency 
range v ft S v S v b is 

^ (x, v)dv 



«. b (X) 


(2.J-L5) 
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where q^U, v) is given by equation (1.32). Using equation (2.112) in 
equations ( 1 . 26 ) and (1.32) aad reversing the order of integration in 
equation ( 1 . 32 ), equation (2.113) may be written as 


Q ab (x) -2 n j sgn (r(x) - t(x' )U ab (x,x' ) dr'(x') 


(2.114) 


where J (x, x* ) is the following integral 
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( 2 . 115 ) 


and t(x) is an optical length which is independent of frequency. 


. . / dx' 

t(.x) = j 

^ O(x') 
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( 2 . 116 ) 


Equation (2.114) for the radiative heat flux is in the form expected 
for a grey gas where J is to be interpreted as a combined attenuation 
function and source function. 

Consider the frequency range £ v < ». Noting that 

hv/KT e « i, equations (2.110) and (2.111) give the following relations 

6<x > * ;cV 

a o 


^ - KTJ 

1 - CUT ) 2hv 5 -hv/4T (x) 

S(*,v) - -- i-a — — e 


Using these formulas the function J(x, x' ) may be written as 
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Jljx, x' ) • v x YZ 


- Og(T e ) 


B (T ft )J du JdA(HA) 


-bVj A - Lt(x) - r(x»)J > 
, kT e u( 1+A) 2 


where A is defined by 

V = V 1 (i + A) 

Since ( hv^ / kT g ) is large, the main contribution to the 
integral will be for small values of A. Assuming A « 1 and 

hv l 

pjr » |t(x) - t(x’)| , J^(x, x ') simplifies to 

e 


J (h . X 1 ) = T" 

-Loo 9 7 1 - a 


1 - a (T ) kT 

- ”iT^(T e ) e g (| T (x) - t(x')|) 


liquation (2.114) then gives the radiative heat flux in the frequency 
range < v < » as 


kT 


, 1 - a.(T ) 

= 211 J J a ‘ — — (Tj IT 3gn " T ( x ' )l E 2 r - |t(x) 


. r(x' ) | ) dr(x' ) 


(2.117) 


This expression is in the form of tne radiative heat flux lor a 
grey gas with an equivalent grey source function 


1 - Cl.(TJ 


a' e‘ 


1 - Q 

and optical thickness 


kT 

B (T ) r-— 

e' h 


x 

I dx • 


T(x) * J ITRTb, 


(2.11b) 


A similar reduction can be carried out for the frequency range 

v u < v < v,. Assuming the Gaunt factor G and effective charge Z 
M — jL 61 1 
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to be one, and noting that — — » i, the following relations are 

K e 

found from equations (3.110) and (2.111). 


B (T e ) JT 

M e 


*-)- 

- v m/ 


3 -hv/nT 

e/ „ V 2ny ^ ' e 

&(*,v) a ~T e 


Substituting these relations into equation (2.114) gives the radiative 
heat flux from the frequency range v M < v < as 

* kT 

X (x) = 2n j B v (T 0 ) ~ agn It'(x) - t' (x* ) Jfigl |t' (x) - t' (x* ) Jd-r’ (x' ) 

■°° (2.119) 


where t'(x) is the following optical thickness 


I K ' n e (x '> 

T*( X ) » / T dX' 

J o B y [T e (x')l T*(x') 

M 


( 2 . 120 ) 


For frequencies between 0 and v., the penetration length is 

ri 

not exactly separable into a function of x times a function of v because 
the combination occurs in the PIsrck function. However, the 

penetration length may be written in an approximate form which i6 
separable. The penetration length in this frequency range is 


t (x,v) ** 


• 4 T B (T ) 

e v e 

- hv/kT 

K n 1 - e 

e 

^ T e B M (T e } 

: - hv/kT 


K' n. 


1 - e 
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where B^(T e ) is the maximum value of the Planck function at the 
electron temperature T g . It is a function only of the electron tem- 
perature and is given by 


£(kT e ) 5 

B M (T ) - 1.U2H06 

M e (he) 2 

The ratio of the Planck function to its maximum value is 

tiv 

which is a function only of r— = — . 

K i e 


( 2 . 121 ) 



/hv A V T e> 

r *v * y 7 ? 


\ f* 

\ 

/ 

1.42145b hv/kT 

e - 1 



Approximately, the penetration length is 


J(x,v) 



\ <V 


K’ 


n~ 

e 


1 



- e 


h v \ 
K>r av J 
-hv/k T a; 


( 2 . 122 ) 


(2.125) 


where T av is a suitaoly chosen average temperature. 

Using the approximate penetration length the following relations 

are obtained. 


e(x) 


4>(v) 


K w 
K ' 4 

r( hv/ hT ay ) 

^hV/kT " 
, ' av 

1 - e 


S(x,v) = B v (T e ) = B fl (T e ) r ( 


hy 

kT 


av 


J 


Substituting these formulas into equation (2.115) gives 
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o o 


kT 

- M iTe (^) 1 F[|t"(x) - t m (x»)|] 

hv M 

where, noting that — - = T gx , the function F may be written 

T /T 
ex' av 

F(t) * j d£(l - f t 

o 

and t" is the following frequency Independent optical thickness 

, X K' n 2 (x’) dx' 

T " ( X ) = j h (2.125) 

o ^T e (x* ) B M [T e (x')J 

The contribution to the radiative heat flux from the frequency 
range 0 i v * v M is therefore 

00 kT 

* £ " j W iT^ w [t “ (x) • T~(x*)iF[ |x~(x) - t”(x')|] dV'(x') 

(S.126) 


,1 - e 

nrr 


(2.124) 


Although equation (2.126) lias the form of a radiative heat 
flux for a grey gas, it does not have the attenuation function 
appropriate to a grey gas. Instead the radiation is attenuated 
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according to the function F. 

The attenuation function F wae evaluated on an IBM 7094 computer. 
The results for T fty equal to ten thousand and thirty thousand degrees 
Kelvin are shown oa Figure 2.1. The grey gas attenuation function is 

also shown for comparison. F is much more sharply peaked for small 
values of its argument than is . Both functions rapidly diminish as 
their arguments increase. At the larger values of the arguments shown 
in the figure, F becomes independent of T qv . 

The differences between che functions F and show that the 
widely practiced procedure of using frequency averaged penetration 
lengths in grey gas radiation energy transfer equations in order to 
approximate the radiation energy transfer in a non-grey gas may be 
erroneous because the grey gas attenuation function may t»e improper. 

The total radiative heat flux is given by the sum of equations 
(2.117), ( 2 .H 9 ), and (2.126). Dividing the total radiative heat flux 
by p Q u^, which i6 the way in which it appears in equation (2.66), 
and expressing the result in terms of the non-dimensional variables, 
the following formulas for the radiative heat flux are obtained. 




3 

P « 
o o 


where T = 


r 1 % M (x) + Q Mi (x) + QjJx)] = r q K 


( 2 . 127 ) 


0 u3 
M o 0 


is a characteristic ratio of radiated energy to 
50 


convected energy and y* ~ (1.421436)> y is the Stefan- 

n 

Boltzmann constant* The separate terms in this equation are 
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, / T \ 4 

Q li( (x) m jg J Pi 1 " C» B (T # )]r ^ ° n - j T e 8gu[T(x)- T (x' ))B g l|T(i) 

— T e - T (x' ) I ]dx' (2.126) 


5^(5) - J~ J a^.Jr^gn (t'(x)- ▼'(x^lBgllT^xJ-T'tx'JIJdx' (2.129) 


5 oM (x) - j- j a S ?jf e 8gn[T"(I)-r*(x') m|r"(x)-r"(i>)|)<lx* (2.130) 


In these equations l is a characteristic length for the ground state 

& 

photoionization process and is a characteristic length for photo- 

ionization from excited states. 


I = = — 

8 — B 

m o 

a 


£. 3 
eA 


/t o V T o> 


K'(0> a )‘ 


( 2 . 151 ) 


In terms of the nondimensional variables the optical thicknesses 


are 


II 

i- 

i a 

3* 

i 

H 

« o 


T '(i; * 1 j 

ex J 

* 2-2 
1 2-4- 

o rfc 

e 


x 0 „ 

T "W * f- j 

ex J 

l Q - p , 

?f 7/2 


dx 


r (T ex /T e ) 


( 2 . 132 ) 


(2.133) 


(2.134) 

7* 

o 'e 

It should be noticed that these equations depend only on the 


electron temperature, not the atom temperature. Thermal nonequilibrium 
of the gas will therefore affect the radiative flux. 
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The Radiative Ionization Rate 

The rate of radiative energy absorption per unit volume of gas 

4Q»(*,v) 

In the frequency range v, v ♦ dv is - ^ — . Since the line 

radiation naa been ignored, all the radiative energy that is absorbed must 
produce ionizations. Each absorption process, that is, each ionization, 
requires energy hv. The total rate of production of photo-electrons 
is therefore 

a - y*.v> 

* • ST J ~ — av 

o 

*Lxccpt for the factor , tne integral in tnis equation is Just like 

that in equation ( 2 . 115 ). Consequently the procedures used in the last 
section to evaluate the integral in equation (2.115) ere applicable 
this integral. In terms of the nondimensionai variat '.es the radiative 
ionization rate ratioed to the characteristic convective rate N is found 
to be 


" i d %< x > v) 

J ^ — — dv 




L T. T T J 

ion ex av 


( 2 . 135 ) 


where ^ j is the same as 4i ^ given in equation (2.150) except that 
the attenuation function F is to be replaced by the function F* written 


below. 


00 



O 


r"(x')|]dx' 

( 2 . 136 ) 

( 2 . 137 ) 
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Numerical evaluation of F' shows it to be similar to the function 
F as may be seen in Figure 2.2. 

Method of Solution 

The radiation flux and the radiative ionization rate are in the 
form of integrals which have integrands containing variables which are 
not known until the problem has been solved. Yet the problem may not be 
solved until the radiation flux and radiative ionization rate are known. 
This uiiemira is avoided by using u starting solution to calculate tne 
radiative terms. The equations are then solved using the radiative 
terms determined by the starting solution. The resulting solution is 
useu to re-evaluate the radiative terms which are then used to generate 
still another solution and so forth. 

The starting soluticr. which was used is that obtained by omitting 
all radiative terms from the equations. The radiative terrs determined 
from this radiationless solution may not be directly used in the equa- 
tions to generate another solution because this procedure does not con- 
verge to the true solution rapidly enough to be practical, instead the 
radiative terms obtained from the radiationless solution are inserted 
in the equations and multiplied by factors which make them artificially 
small but not negligible. The factors are chosen so that the next 
solution does not appreciably differ from the starting solution. This 
procedure lb repeated, each time increasing tne factors which multiply 
the radiative terms. In this way a series of solutions are generated, 
each of which is not greatly different from the preceding solution. 

When the factors are one the true solution is obtained by iteration. 
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Knowing the values of the variables at an initial point x^, 
the values of the variables at any point x may be found as follows. 

Using equation (2.133) the electron mass conservation equation 
(2.69) becomes 


, r ftA > * ft (i) 

Sti ( i> - ftft 

fti ft) - 1 

" a i + a - + 

4 T. 

ion 

T 

ex 

J 

av 

X 



*J.( — ) a* 

x, N / 


( 2 . 158 ) 


where 0 ^ is the initial value of the degree of ionization and the terms 
StL> anci ^oM are to be re 8 arde( i as Known functions of x.from the 
preceding solution. The integral in equation (2.158) was evaluated 
using a standard Runge-Kutta integration procedure. The integrand of 
this integral is found from equations (2.85) and (2.84). It must be 
evaluated at the values of a and x generated by the Runge-Kutta 
method. This is accomplished in the following way. 

Equations (2.64 ) and (2.6o)may be combined to give 


p(x) = 


5C. 


'“ft ’ ° 5 io»- l % 1 


I 

L 


1 + 


(x- 


J2.(c e - cJ iOB - £ q„hi 1 


250' 


-TJ 


- = (c 

T + a r P 

a 6 


i) -ft 

p a p 


(2.159) 
( 2 . l40 ) 


The density p can be calculated from equation (2.159) using the 
values of a arxl x generated from the Runge-Kutta integration proce- 
dure (q n is a known function of x from the preceding solution ). 

Then equation (2.140) can be used with the electron energy equation 


The average energies ^ 8111(1 T rad which appear in the electron 
energy equation are shown in Appendix C to be of the order kT or less. 

(B 

When the electron-atom rate A^ is the dominant rate, the terms la the 
electron energy equation which conte in the average energies are com- 
pletely negligible. Just behind the viscous shock front the atom-atom 
rate will be the dominant rate. There will then be some error in 
determining the electron temperature if the terms containing and 

^rad are omlt ' be<1 ^ rom the electron energy equation. However, when the 
atom-atom rate is the dominant ionization rate an accurate determination 
of the electron temperature is unimportant because then the atom-atom 
ionization rate is independent of the electron temperature. In the 
precursor region, where only the radiative rate is important, the 
electron temperature need not be calculated because the degree of ioniza- 
tion in the precursor region is so small that the terms containing the 
electron temperature are eliminated from the equations. Consequently 
the terms containing the average energies and £ ^ were omitted 

from the electron energy equarion. 

The constants Cp and Cg are^ of cour6q determined from the 
initial conditions at point x^. 

Theoretically the initial conditions are known at x ■ - «o ahead 
of the shock wave. To start a numerical solution of the problem, however, 
initial conditions must be known at a finite distance in front of the 
shock wave. In order to find the conditions at a finite distance ahead 
of the shock wave from those which are known at - <», an analytical 
solution valid far ahead of the shock wave is used. 
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In the precursor region far ahead of the shock vave the gas will 
be cold. The temperature will be so low that the radiative ionixation 
rate will be much larger than the colllsional ionixation rate. Further- 
more, the number of excited atoms will be small and consequently the 
radiation emitted by electron recombination to excited states in the 
shocked gas will not be absorbed. Accordingly equation ( 2 . 69 ) can be 
integrated to read 
00 

m a I (Jy 

a + ■ I — q H (x,v) = constant (2.141) 

p u J hv 

K 0 O 


In the far precursor region the emitted radiation will be negligi- 
ble compared to the absorbed radiation. Consequently equation (1.24), 
the radiative transfer equation, becomes 


» ar ixl v t~ v J - - - I(I > v ’ “> (2 - iW) 

The radiation intensity at x = 0 is found from equation (1.57) to be 


1(0, K-,v) 



S(x* , 


V ) e -T (x,v)/M _dT 'x',v) 
d 


if -l£ d £ 0 
if 0 < d £ 1 

(2.145) 


The solution of equation (2.142) with the requirement that the solution 
at x = 0 be given by equation (2.145) is 

I(x, d, v) = 1(0, d, v)e ' t(x ' v)/u (2.144) 

Using the Eddington approximation the radiative flux can be 
related to the average intensity as follows. 
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,v) ■ 2« J (i I(x, u,v) d|i » 4* Hi (x,v) 

-1 


(2.145) 


where 4 Is an average value for \i, assumed to be and the average 

intensity, I, is defined by 

1 

— 1 [ 

I(x,v) = -g J I(x, U, v)d U (2.146) 

-1 

Equation (2.142) may be integrated over to give 


dq K* x *v) _ , dl( x >v) 

dr(x,v) ” n dr(x,vT 


4n I(x,v) 


The solution to this equation is 


I(x,v) = I ( 0, v) e ' T ( x > v )/^ 


(2.147) 


Collecting these results the integral in equation (2.141) may be 
written as 

-t(x,v)/u f 1 


m t * m , . -t(x,v)/v r 

P Q « 0 J *V ^R (x > v) = p o U Q hv J 


du . 


-t(x»,v)/u d / , ) 

.J S(x',v) e t f" ^ 


o 


Using the penetration lengths and source function appropriate to fre- 
quencies greater than from equations (2.110) and (2.111) this may be 
written in the form 


m 


P u 
H o o 


J \ 


(2.146) 
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In this equation all quantities have been referred to the knovn initial 
conditions at x » - « 


P 0 - P( — ) # - T( -^), u Q • u(— ), etc. 

except T and the variables in the integrand of Q(o) which are referred 
to values typical of the hot shocked gas downstream from the precursor 

region. Using I for the characteristic length J, the optical thickness 

& 

which appears in equation (2.148) may be written as 


T(X) 


J 


(1 - a) p dx 


(2.14 i. 1 ) 


In a similar way, the radiative flux is found to be 


-t(S)/p r - 

^(x) = 4n j I (x,v) dv - e J. (o) 

V 1 


The degree of ionization in the precursor region far ahead of 
the shock will be negligible compared to one. Since the electron 
temperature always appears in the equations multiplied by the degree 
of ionization, the electron temperature may be omitted from considera- 
tion in the far precursor region. The equations become 


- + a p T = C 

- M a p 
P 

*5 — — i r ■ 

^ T_ + a T . _ 4 + — e 


t/u 


2 a 


i0n 2a r a 


QlJ°) - C g 


p “t / d Hi \o) 

a + — e ^ — a constant 

T 

ion 


At x = - «, p = T =1. The solution to these equations is 

& 


(2.150) 

(2.151) 

(2.152) 


4 
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p(x) « T a (x) - 1 ( 2 . 153 ) 

-t(x)/H 

a - a(o) e (2.154) 


a(o) 


a T 


ion 



( 2 . 155 ) 


This solution gives the finite value of x needed to start the 
solution and the initial values for p, T and a in terms of the 
ratio a/a(o) at which the solution is to be started and the flux 
integral 5 ^( 0 ). 

The full equations are then integrated, using the procedure 
previously discussed, to the point x = 0, where the gas passes through 
a shock wave caused by viscosity and heat conduction effects. The 
radiative flux is constant and no ionization takes place across the 
viscous shock wave. The fluid conservation equations then give the 
jump conditions. 

p(o + ) = 4 0 ( 0 ") 


q(o + ) = a(o”) 


T a (o + ) + o T e (o + ) = 


“ u*(- -) 


p(o“) 


1 

t 


where the - and + signs denote values Just before and Just after 

the shock respectively. The la6t equation must be solved with the 

, + . + 

electron energy equation to get T (o ) and T (o ). 

a e 

From behind the shock the equations are integrated until 


equilibrium is reacied to within 6ome predetermined accuracy. The 
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definition of equilibrium used for the calculations was that the degree 
of ionization equal the local equilibrium degree of ionization within 
five percent and the electron and atom temperatures be equal within one 
percent. The equations are very difficult to integrate near equixibriun 
because then the ionization rates are sensitive to small changes in the 
variables. 

After equilibrium has been reached the degree of ionization can 
be calculated from the Saha equation, and the system of equations 
becomes entirely algebraic. The variables change in this region because 
of radiation cooling. The radiative cooling continues until a steady 
state is reached. 

In this way a new solution is generated from the preceding solu- 
tion. Tlus process is repeated, each time increasing the factor which 
multiplies the radiative terms, until the factor is one and the solutions 
converge. Successive steps for the Mach number It) solution presented in 
Chapter III are shown in Appendix D. 
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CHAPTER III 

RESULTS AND DISCUSSION 

Computer codes using the method of solution outlined in the 
previous chapter were developed and calculations were carried out on 
on IBM 709 U computer. 

The argon gas ahead of the shock wave was taken to be at a 
pressure of 1 cm Hg and '}00*K temperature. Calculations were performed 
for shock waves having Mach numbers of 12, 18, and 30, which correspond 
to shock velocities of 3 . 87 , 3.61 and 9.68 mm/p sec respectively. 

The reference state used for the calculations was that given by 
the strong shock Rankin® -Hugoniot relations applied to the cold gas far 
ahead of the shock wave. The reference values found by this procedure 
ere shown in Table 3.1. 

Table 3.1. Reference Conditions. 



p 0 U° 5 ) 

(gm/enr) 

u 

0 

(%/n eec) 

T e (l 0 ‘ 4 ) 

# K 

Q o 

r 

a 

M = 12 

1.090 

.7581 

I .096 

.0344 

7.555 

5.968 

M = 18 

17.74 

.6584 

1.423 

.2028 

16.91 

6.076 

M = 50 

25.29 

.8167 

1.987 

.7487 

28.11 

6.202 


These reference conditions give the characteristic lengths shown 
in Table 3.2. 
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Table 3.2. Characteristic Lengths 



1 AA 

(cm) 

(cm) 

V“ 2) 

(cm) 

'ex 

(cm) 

M * 12 

2.387 

1.505 

1.759 

1.791 

M a 16 

.9112 

.498L 

1.068 

1.693 

M = 30 

A536 

.2480 

.7*9* 

2,67* 


At each Mach number solutions were first obtained by neglecting 
thermal nonequilibrium and radiation effects. Then the calculations were 
again performed with thermal nonequilibrium, or radiation, or both included 
in the analysis. In this way the effects of radiation and thermal non- 
equilibrium can be evaluated by comparing the solutions obtained when 
radiation and thermal nonequi librium have been included with the solutions 
obtained when radiation and thermal nonequi librium have been excluded. 

Radiationless. Thermal Equilibrium Solutions 

When radiation is omitted from the analysis and the electron and 
atom temperatures are assumed to be equal the solutions shown in Figures 
3.1, 3.2 and 3.3 are found. 

The temperature and density are to be r<-ad from the linear scale 
on the left. The electron-atom and atom-atom ionization rates arr. shown 
divided by the total ionization rate Just after the viscous shock, n^. 

The ionization rates are to be read from the logarithmic scale at the 
right as is the ratio a/u o . 

'Two scales are used on the abscissa. For the extent of the first 
scale the atom-atom ionization rate exceeds the electron-atom ionization 
rate and the characteristic length is chosen to be 


For the extent 


82 



Figure J.l. M = 12 Radiationless Shock Wave with T 


2.0 
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of the following scale the electron-atom rate exceeds the atom-aton rate 
and the characteristic length is taken to be f The portion of the 
shock vave where the atom-atom ionization rate is dominant is small 
compared to the portion where the electron -atom rate is dominant. The 
length of the region in which the atom-atom rate is dominant, x^, and 
the length of the following region in which the electron-atom rate is 
dominant, x^, i6 shown in Table 5 . 3 . The region where the e-A process 
is dominant is terminated when the degree of ionization is within 5 ^ of 
the local equilibrium degree of ionization. 

Table 3*3. Length of A-A and e-A Regions 

for Radiationless Shock Wave with T * T . 

e a 



X AA (C “> 

x eA (cB) 

M = 12 

•9571 

36.3 

M = 18 

6 . 38 ( 10 ”^) 

8 . 71 ( 10 " 2 ) 

M = 30 

2 . 72 ( 10 " 5 ) 

1 . 92 ( 10 " 5 ) 


In all cases the e-A ionization rate quickly rises to a peak 
value and then slowly diminishes. The peak and decrease in the rate is 
due to the Increased importance of recombination as an equilibrium degree 
of ionization is approached. The A-A ionization is unimportant except 
right behind the viscous shock. 

The behuvior of the other variables follows that of the e-A 
ionization rate, a region of abrupt change followed by a much slower 
variation. The shock waves all have long tails where the values attained 
by the variables are close to their final values. 
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Thermal Equilibrium with Trapped Kadiation 

For the solutions shown In Figures J.k, 3*5 and 3*6 the only 
radiation process allowed was photoionization from the ground state. 

This radiation Is completely absorbed in the shock wave, i.e., the 
radiation is trapped. As la the last section the gas is assumed to be 
in thermal equilibrium and the electron temperature Is equal to the 
atom temperature . These same assumptions were used in the solutions 
obtained by Clarke and Ferrari [ 5 ). 

Ionization takes place in the precursor region due to the absorp- 
tion of the radiation omitted by the gas behind the viscous shock at 
x = 0 . The degree of ionization in the precursor region is shown divided 
by the degree of ionization at the viscous shock wave, , in Figures 
3 . 4 , 3«5 and 3.6 and is to be read from the logarithmic scale at the 
extreme left. The abscissa in the precursor region is x divided by 
the ground state characteristic length evaluated for conditions at 
x - co, t, - .08878 cm. The precursor degree of ionization is shown on 
a larger scale in Figure 3.7 where it can be seen that the variation of 
the degree of ionization is nearly exponential except close to the 
viscous shock wave and in the far precursor region. 

Because of the precursor ionization the atom-atom rate behind 
the viscous shock at x = 0 is much less than the electron-atom rate. 

In contrast to the radiationless solutions discussed in the last section 
there is no region where the atom-atom ionization rate is the dominant 
rate. The peak value of the electron-atom rate is reduced and slightly 
shifted, but it 6tlll lias the same features as for the radiationless 
case, a rapid rise to a peak value followed by a region where it 
slowly decreases. 






5 


-X/lg. 


The distance required for the degree of Ionization to come within 
of its local equilibrium value is only slightly different than when 
trapped radiation Is excluded from consideration as In the last section. 
The distances are 56 . 3 , 8.97(10”^) and 1.84(10"^) cm including trapped 
radiation compared to 57*3, 8.77(10* ) and 1.95(10”^) cm when trapped 
radiation is excluded, for Mach numbers 12, 18 and 30 respectively. 

Thermal Nonequilibrium with Trapped Radiation 

As in the last section the only radiation process which is con- 
sidered in the calculations is photoionization from the ground state. 
However, the gas is not assumed to be in thermal equilibrium and the 
electron temperature may be different than the atom temperature. The 
solutions obtained by employing this model are shown in Figures 3.8, 

3.9 and 3.10. 

The electron temperature is much lees than the atom temperature 
right behind the viscous shock. It then rises and becomes equal to the 
atom temperature farther downstream. 

The electron-atom ionization rate depends principally on the 
electron temperature. The lower electron temperature causes the 
electron-atom rate to be considerably reduced from the value it would 
have if the electron temperature were equal to the atom temperature. 
Because of the lower e-«A ionization rate, the point at which the e-A rate 
attains its peak value is much farther downstream from the corresponding 
point when thermal equilibrium is assumed. Also the length of the 
region where the atom-atom rate is dominant, x^, and the length of 
the following region where the electron-atom rate is dominant, 
are ouch greater than the corresponding lengths when thermal equilibrium 









Trapped Radiation Shock Wave vith 
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is assumed. These lengths are shown in Table 3.4. 

Table 3.4. Length of A -A and e-A regions for 

Trapped Radiation .">hock Waves with T + T . 

© © 





x eA (cm) 


M « 1? 

7.17 

b4.2 


M » 10 

.?05 

• 3**9 


M - 30 

3.lb ( 10 ' 2 ', 

3.47 U0~ 2 ) 


Because of the lower electron temperature the intensity of the 
radiation emitted fron the gas just behind the viscous shock will be 
significantly iess than the intensity of radiation which would be .-mitte.i 
were the gas in thermal equilibrium. The penetration length for ground 
state continuum radiation la less than or of the order of the length re- 
quired for the electron temperature to become equal to the atom tempere- 
,uxe and all of the ground state continuum radiation is omitted into the 
precursor region from right behind the viscous shock. Consequently tne 
Intensity of ground state continuum radiation is reduced. The precursor 
ionization for Mach number 12 is so much curtailed it is negligible. For 
the Mach numbers 10 and 30 reduced precursor Ionization occurs. The 
degree of ionization produced by the Mach numbers 18 and 30 thermal non- 
equilibrium shock wave is shown in Figure 3.7. 

The preceding results show the ratio of the penetration length to the 
length required for the electron temperature to attain the atom tempera- 
ture is an important parameter in analyzing the precursor ionization 
caused by emission of ground state continuum radiation. It can also 
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be an Important, parameter in determining precursor ionization due to 
emission of radiation caused by other radiative reactions involving 
electrons. 

Since the calculations show precursor ionization due to emission 
of ground state continuum radiation is negligible at Mach number 12, and 
considerable precursor ionization is measured at Mach numbers and 
pressures similar to those used in the calculations, [37], (381, [40], 
precursor ionization must be caused by some other radiative reaction. 

Thermal Nonequiiibrium with Complete Radiation 

The analysis of this section includes radiation resulting from 
photoionization and radiative recombination involving excited bound states 
of the atom in addition to the radiation due to photoionization and ra- 
diative recombination involving only the atomic ground state which was 
considered in the previous section. The equations then include the 
integrals specified by equations (2.129) aad (2.130). 

The gas is not assumed to be in thermal equilibrium and the 
electron and atom temperatures may be different. Since the radiation 
integrals in equations (2.129) and (2.130) depend on the electron 
temperature, as do the equivalent optical thicknesses given by equations 

(2.133) and (2.134), the contribution of these terms is affected by 
the thermal nonequilibrium of the gas. Equations (1.26), (2.133) and 

(2.134) show the effect of lowering the electron temperature, all other 
factors being the same, is to decrease the penetration lengths of the 
excited 6tate continuum radiation. Equations (2.129) and (2.130) show 
the excited state continuum radiative flux will be reduced if the 
electron temperature is lowered. Accordingly the radiative flux and 
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penetration lengths of the excited states continuum radiation is changed 
by the cool layer of electron gas Just behind the shock discontinuity 
when thermal equilibrium is not established. 

The solution for a Mach number 16 shock wave is shown in Figure 
5.11. It consists of an atom-atom ionization region and an electron- 
atom ionization region, just like the corresponding solution without the 
excited state continuum radiation shown in Figure 5*9» followed by a ra- 
diation cooling region. The length 1" chosen to aondimensionalize 
the coordinate x in the radiation cooling region is a typical penetra- 
tion length for the continuum radiation in the frequency range 
0 § v £ v .. It is the distance in the uniform slab which would result 
after the shock if there were no radiation cooling for the optical 
thickness r" to be one. For the Mach number 13 shock wave 1" * U2.7 
cm. The optical thickness t' corresponding to continuum radiation in 
the frequency range v u S v S v, is always much larger than r" and 
the radiative flux term is negligible compared to 

Up to the beginning of the radiation cooling region the solu- 
tion is almost identical to that found excluding the excited states 
continuum radiation. At the beginning of the radiation cooling region 
the solutions for these two cases are within l^t of agreement. The 
equilibrium degree of ionization, however, is very sensitive to small 
changes in the solution and at tne beginning of the radiation cooling 
zone, where without the excited states continuum radiation the solu- 
tion was within of the equilibrium degree of ionization, the degree 
of ionization is only within about y)f> of the equilibrium degree of 
ionization. Furthermore, 86 equilibrium is approached, the collisional 




lonlzatlonal rate becomes small and the radiative rate becomes relative- 
ly more Important. The radiative rate Is a recombination rate. It 
opposes the colUslonal ionization rate. As a result the point at 
which the degree of ionization reaches 5 ^ of its equilibrium value is 
considerably delayed. For the solution shown in Figure 3.11 this point 
is at x/i" ■ .1045 which is 4.46 cm from the end of the e-A region 
where the equilibrium point was attained when excited state con- 
tinuum radiation was excluded from the analysis. The end of the e-A 
region i 6 .345 cm. from the shock wave discontinuity and corresponds 
to 5.38 usee. laboratory time. This compares well with the time to 
reach equilibrium ionization measured by Petschek and Byron [ 19 ) of 
6.5 usee, and that measured by Wong [20] of 4.7 usee. The analytical 
solution of this section suggests, however, there is an extensive region 
beyond this point where the gas has not attained equilibrium, but is 
near equilibrium. 

The distance to reach an equilibrium degree of ionization cal- 
culated in this section is that corresponding to T = 30,000 *K. This 
is a conservative estimate. Smaller values of T will result in longer 

ft V 

distances to reach an equilibrium degree of ionization. Since line 
radiation, which has been excluded from tne analysis, can cause an 
increased ionization rate and shorten the distance required to reach 
an equilibrium degree of ionization, the extent of the chemical non- 
equilibrium region is uncertain. It is clear, however, the relative 
magnitudes of colJLislonal and radiative ionization rates are important 
when chemical equilibrium is approached. 
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CHAPTER IV 

SUMMARY AND CONCLUSIONS 


A gas, initially composed only of atoms, becomes a mixture of 
atoms, ions, and electrons in the process of ionization by a shock wave. 
The ways in which mass, momentum and energy may be transferred within 
such a mixture are different than for the unionized gas. The shock wave 
is profoundly affected if the transfer processes occur slowly compared to 
ionization. This is indeed the situation for the shock waves investigated 
in this thesis. The electrons lose translational energy in the process of 
ionizing neutral atoms. Furthermore, since electrons are inefficient in 
exchanging energy with the atoms and ions, they regain energy slowly, 
being effectively insulated from the atoms and ions. On the other hand 
the atoms and ions readily exchange energy. In this situation the gas 
may be thought of as a combination of two component gases, a gas composed 
of only electrons and a gas composed of both atoms and ions. The tempera- 
ture of the electron gas may be different from i.hat of the atom-ion gas, 
that is, the electron gas and the atom-ion gas may not be in thermal 
equilibrium. The inability of the electron gas and the atom-ion gas to 
establish equilibrium greatly affect6 the rate of ionization and radia- 
tive emission behind the shock wave. 

Viscosity and heat conduction are changed by ionization. The 
electron gas is an especially good conductor of heat. For the shock 
waves of interest in this thesis, viscosity and heat conduction effects 
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are omly Important in a very thin portion of the shock wave. In this 
case the effect of viscosity and heat conduction can be thought of as 
dlscontinuously changing the temperatures and density at a point in the 
gas. Following the discontinuity the atom-ion gas cools because it 
loses energy by creating further ions, heating the electron gas and by 
emitting radiation. The temperature of the electron gas, however, is 
less than that of the atom-ion gas after the discontinuity. Subse- 
quently there is a net gain of energy by che electron ga3, heating it 
until thermal equilibrium is established with the atom-ion gas. 

Two collisional processes are responsible for ionization, the 
atom-atora process and the electron-atom process. The rate of ionization 
due to the atom-atom process depends on the properties of the atom-ion 
gas, but the rate of ionization due to the electron-atom process depends 
on the properties of the electron gas. Since the electron gas is 
cooler than the atom-ion gas, the ionization rate due to the electron- 
atom process will be retarded compared to the rate which would be 
predicted assuming the electron gas and the atom-ion gas were Ln thermal 
equilibrium. The results presented in Chapter III show the ionization 
is significantly different in these two cases. 

The grey gas assumption implies that radiation energy transfer 
can be described by a single frequency independent penetration length. 
For argon, however, there are many different penetration lengths because 
there are many possible radiative interactions which can occur. Conse- 
quently a grey ga& model for argon is open to question. 

The excited states of argon are closely spaced within a few 
electron volts of the ground state ionization limit. It follows that 
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the energy levels may be approximated as continuously distributed and 
hydrogenlc. The radiation energy transfer due to photoionization and 
radiative recombination of all the excited states may then be found by 
integrating over the continuously distributed levels. The resulting 
radiation transfer equation has a marked frequency dependence. It 
agrees with the Unsflld result for small values of frequency, but at 
large values of frequency, where the Unsold model fails, it has the 
proper form. For the shock waves considered in this thesis the pene- 
tration length for photojonization and radiative recombination from the 
excited states is large compared to lengths typical of other processes 
which may take place. In the precursor region it may be regarded as 
infinitely large. Consequently continuum radiation caused by photo- 
ionization and radiative recombination of the excited states results 
only in the radiation coolin tJ of the gas behind the discontinuity. The 
penetration length for photoionization and radiative recombination of 
the ground 6tate, however, is sufficiently small that the continuum 
radiation from the ground state is trapped, causing precursor ioniza- 
tion. 

The total radiative flux and the radiative ionization rate are 
in the form of multiple integrals. These Integrals can be reduced to 
the form they would have if the gas were grey, but the source functions 
and attenuation functions contained in these integrals in general are not 
those of a grey gas. Consequently the practice of using grey radiation 
transfer formulas to describe radiation energy transfer in a non-grey 
gas is not Justified for the shock waves investigated in this thesis. 
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Mediation energy transfer depends on the behavior of the electron 
gas. Just after the shock discontinuity the radiative flux is less than 
the value it would have if thermal equilibrium were established, because 
the electron temperature there is less than the atom temperature. For 
the calculations performed in this thesi6 all of the radiation which is 
trapped in tne precursor region comes from the cool layer of electron 
gab just behind the shock discontinuity. Consequently the precursor 
ionization \g much reduced from that, which would be estimated if thermal 
equilibrium were established. In fact the precursor ionization for the 
Much number 12 shock wave is found to be negligioly small because thermal 
equilibrium is not ectabiionea. Consequently the measured precursor 
ionization [37], [ 38 ], l 40], for Mach numbers near 12. is not caused by 
continuum radiation trapped in the precursor region. 

The penetration length for radiation caused by pnotoxoniza tion 
from the ground state 16 independent of temperature, but the penetration 
length for radiation caused by photoionization from excited states depends 
on the temperature of tne electron » as. The excited state penetration 
length is an equivalent penetration length for photoionization for all 
excited levels and all possible frequencies of radiation. As shown by 
equation (2.111) it depends on the nurnuer density of electrons and is 
proportional to the electron temperature. As a result the penetration length 
for excited state continuum radiation is much smaller in the ccoi layer 
of electron gas near the shock wave discontinuity than it would be if 
the thermal equilibrium were established, fnis layer is therefore much 
more absorbent when thermal equilibrium is not established m the 


shock wave 


lJU 

Radiation energy transfer, therefore, depends on the relative 
extent of the thermal and chemical nonequilibrium regions compared to the 
penetration lengths of radiative processes. The extent of the thermal 
and chemical nonequilibrium regions Is reduced when the Mach number is 
increased. 

Radiation has little effect on the ionization until equilibrium 
ionization is approached. Then the collisional rates decrease and the 
radiative rate becomes relatively more important. Since the radiative 
rate opposes the collisional rates the effect is to delay the point at 
which equilibrium ionization will oe reached. There is an extensive 
region where the gas properties are nearly constant but not in chemical 
equilibrium. Since photoexcitation processes, which could increase the 
ionization rate and cause equilibrium to be reached more quickly, have 
been excluded from the calculations, the existence of the quasi- 
equilibrium region i6 questionable. It i6 certain, however, that the 
rate of ionization near equilibrium ionization is a delicate balance 


of collisional and radiative ionization rates 
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APPENDIX A 


DERIVATION OF COLLISIONAL RATS CONSTANTS 


The number of collisional encounters per unit volume and time 
between particles of species A having velocities in the range 
v^, + dv^and particles of species X having velocities in the 

range v , v^ t + clv^, is given by the following classical kinetic theory 

expression [4?1. 


n A n A' f A f X' gb db dd dv x dv A' 
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where f is the velocity distribution function 
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b is the impact parameter, 9 is the angle between the plane of motion 
and a reference plane, and g is the relative velocity of the particles. 

ut Hv ‘ (b,g) be the probability that the A particle is 
excited from state 1 to state ? during the collision with the particle. 
The ntm’rer of excitations per unit volume and time is therefore 


Vv Vv z «. 


The to tax numoer of e -citations for all orientations and magni- 


tudes of the vectors v and v* is 
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where S^,(g) is the diffusion cross section 
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end & , s i if species X is different than species X', zero otherwise. 

Consider the collisions between atoms. iixcept for the small 
number of particles that become excited, translational energy is 
conserved in a collision so tha& [42j 
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dv dv , = d& dr 

a a’ ^ 

where G i& clie center ol mass velocity. Substitution it o equation (A. 4) 

gives the number of excitations per unit volume and time as n K.(l,2) 

a * 

where 

r m a n 3/2 r“ p -4 m a 8 2 /itT a 

K A (1 > 2) = 5 ( ?T W J J 4ng ag e g S aa^ ; ’ (A * 5) 

*o 

r- 

is defined by requiring £ n.^ to be the threshold t-nc-rgy 

2 

required for the excitation. Since kT « E < r , the cross 

a ex ■” c 8 

section need only be known near the excitation threshold to obtain an 
accurate value of the integral in equation (A. 5). Substituting 


S 

aa 
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in Equation (A. 5) gives t-.he rate constant used in equation (2.75) 
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k a u,2) *- u v (■ . 8 ) . 
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where T * E /k. The value of C. * 7. 5(10*") cm^ erg used in this 


ex ex 1 


thesis is that determined by Kelley [23] . 

For collisions between electrons and atoms the velocity distribu- 
tion functions are given by equation (A.?) at the temperatures T p and T fi 
respectively. Since the electron and atom temperatures are of the 
same order of magnitude and the elect i’On mass is much smaller than the 
atom mass, the electron velocity must be much larger than the atom 
velocity. This is the basis of the following approximations. 


£ = gP = [v~ - 2v p . v a + v^ p 
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’Jsing these approximations equation (A.J) becomes 


n n / Un v 1 " dv f (v ) S (v ) 
ea/ e ee'e ea e' 


(A. 9) 


where <, m v is the tnreshoLd energy E required for the excita- 


tion. As in the previous result, the cross section need only be known 


near threshold to give an accurate value for the integral in equation 
(A. 9). Near the threshold energy tne cross section may be represented by 
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3 (v ) « C (i m - E ). (A. 10) 

ea' e' e' I 2 e e ex' 

With this approximation for the cross section equation (A. 9) 

gives the number of electron-atom excitations per unit volume and time 

as K (1,2) n n where 
c * e a 

C e 5/2/ V \ - T ex / T e 

K.(l,2) = -2 (2kT e ) ♦ 2) e “ (*.11) 

I a ai 

\ r-» 

\ 

The value of the constant in the excitation cross section formula 
employed by Petschek and Eyron llS>] is C g = 4.4(10 ) cm‘"/erg. This 

value is used for the computations In thi6 thesis. 
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APPENDIX E 

RADIATIVE REACTION TRANSFER EQUATIONS 


It follows from equation (A. 9 ) in Appendix A that the number of 
collisions between electrons with energy in the range q, q + dq and 
ions per unit volume and time is ^ f g (q)dq where f e (q) is the 
energy distribution function for the electrons. 
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I«t A do) be the probability *,hat in such an encounter the 

electron will lose energy hv = X + q oy emittin e ; radiation within a 

P 

A 


solid angle dn) of a direction specified by a unit vector L. The 

th 

electx’on will then be in the p atomic state. The total radiation 
intensity emitted in tue frequency interval v, v + dv per unit volume 
and time by all such collisions is 


hv A n n f (q) h dv dn) (B.l) 

cp ex e 

In addition to the spontaneous emission Just discussed emission 

of radiation may be induced by the radiation field. Let 3 be the 

cp 

probability that in an ion-electron encounter the radiation field causes 
the recombination. The radiation intensity emitted at point P within 

A 

solid angle du) of the direction specified by L and in the frequency 
range v, v + dv induced in electron-ion encounters is 


U5 


hv B c p Lf “e n i f e ^ h dv 


(B.2) 


Let B be the probability per unit intensity that an atom in 
pc 

state p will be photoionized creating an electron with energy in the 
range t[ , rj + dij and causing a loss of radiation energy equal tc hv ■ 

X + Ti from the radiation within solid angle do> about the direction 

p 

A 

specified by the vector L. The total intensity of radiation absorbed 
from the radiation field within a solid angle via) of the direction 

A 

specified by L and the frequency range v, v + dv due to photoionization 
from atomic state p at point P is 


hv n(p) B I(P, L, v) do dv 
pc 


(B.5) 


The rate of change of intensity at point P within solid angle 

A 

da) about the direction specified by L in the frequency range 
v, v + dv i6 the sum of the intensities given by equations (B.l), 
(B.2) and (B.J). 


. - tn( 
ds 


(p) B pc - n e n a f e (r,) h B cp ] h v I( P , L , v) 


+ hv A n n. f (n)h 
cp e i e' •' 


- n(p) B 1 i - f e (n) n ^ lhv 1 (P, L, v) 

* ~ L n( d) oc - 


+ hv F* Tjfpt f e {r]) a 

pc 


(BA; 


Using the Sinstein relations 


i 
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is£ « i 

B po *f 


A 

-S£ 

B 

cp 


2 hy" 

? 


where g is the degeneracy of the free electron states [60] 


5/2 

4* m 

- 2 -T— V 


£J \ 


and is the degeneracy of the p tn atomic state, equation (3.3) 

can be written in i.he form of equation (2.6H) 


dI e (x >“> v) i . 

"3^ ' -J^(x,v) U(< > B > v ' • V A,v) 1 


(8.5) 


vci ere 


, . - n(p) - hv/kT 

— — , ■ = n(p) B I i - -• — • ' e 

V < x > v ) P° L n( 


i(p) 


( B.6) 


n(p) S 2hv 5 / hv / RT e n ^E\ l 


S(x ' v) * TOf ~ I® 


>(p) 


(B.7) 


in tnese equations n(p') is the population of state p which would exist 

E 

if state p were in equilibrium with the ions and electrons. 


n e n i ? ( 2 * m e K ' r e ) 

^ * h 5„ 


Vs 


-1 /kT 

p' « 


(R.*) 


The radiation energy transfer equation for rauiatlve excitations 

can be derived in a similar way. Let A be the probability l’or 

qp 

spontaneous emission, B that for induced emission and b that for 

’ qp PI 

absorption for a radiative transition between atomic states p and q. 

The radiation transfer equation for pnotoe/.citation reaction between 


states p and q is 


— 21 * '"(P^pq P pq ( v ) X ( p » L > v) + ^8)^ P p q(v) + *(q) ^ Ppq(v) 


= - n(p) B (v) P p q(v)[l(P, L, v) - ^ -f ) - ^ j -* ^ ] < 

P q P' 1 L n(p) 6 q n(p) c q c J 


where p^(v) Xs x ^ ne s ^®P e function. This may be written in the 
form of equation 


dI ( x >^ v) 


TTvJ) v) - V x ’ v) 1 

pq 


wlieie 
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APPENDIX C 

AVERAGE ENERGIES OF CREATED EIl£TRONS 


Assume the primary way in which electrons are created in the 

atom-atom ionization process is by collisions; between electrons and 

excited argon atoms. An electron of energy r| on collision with nn 

e /.cited atom will produce an electron of energy £ and lose energy 

K, - X, In the process. 

*: 

e + A(t) ■* A + e - e. y ( C . 1 ) 

*i ' ' n (n - ^ - 0 

where the energies of the electrons are indicated oy their subscripts. 

It follows from Append!.'. B that the rate at which electrons are 
created is 

n a n e J m 1 R 0 dT l d C 

> e 2 c 

where ^ d£ is the cross section for the process described 

by equation (C.l) and f(tj) is the electron energy distribution function. 

-V*T e 


r(i\) 


2n 

(n XTj 






n * 


The total rate of proauction of energy in the form of created 
electrons is therefore 


oo T) - 

V. j fh > ^ d "j * t s 2.c b,t) 15 


1 
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Dividing this by the rate of production of electrons gives the 
average energy of a collisionally created electron. 

C c (n, t) <iC 


J 




J 


X - 




J f(n) T|* QTJ j S 2 ^ c (T|, 


(c.2) 


Using the Thompson cross section 

'Vc (t *> ] = 


ice 


* if 


the integration over £ can easily be performed. The remaining inte- 
gration can be carried out by noting ^ » kT g and expanding the 

integrand about the value t| = X^. 


i 


-n/kT 


[in x: 


- ( 1 - f )] dn 


i - 


J 


-n/kT 


T- ( 1 - 


-V kT e 


“ -n/kT e 

J e 

o 


1 2 ( )* + • • . ] dj) 


, - X„ /’kT 

1 2 ' c 


I 

J 


-n/kT. 


( 3 - + 

1 x„ 


) Urj 


i 


-X^/kT 
c e 


t kT 2 

kT [( ~)+ ... ] 


■X 'kT 


kT e [ ( ^ ) + ... ] 


Consider next the electrons produced by photoionization from the 
ground tate. Assume the radiative intensity is given by the Planck 
function. The number of absorptions producing electrons with energy in 
the range i .» then (c.f. Appendix 5). 


h(l) B. , 

— B (T ) rS 

hv v e h 


The average energy of photoelectrons created by photoionization 
from the ground state is therefore 


i — B (T ) ^ 

s h v v e h 


j r— B (T ) ^ 
J h v v e n 


Making the following, substitution. 


hv - X 


V , ^ 

B, = B ( — ) c 
ic O v ’ 


-nv/kT 

y e 


fC.4) 


B (T ) = ^ 

v e' c 


equation (C.4) can be written as 


i 
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* -hy/kT e 

j (hv - X) e 

V 1 kT 

C — - "3T ( c ‘5) 

I -bv/kT e 
J e dv 


A 3 in Chapter II assume the excited states are closely spaced 
and nydrogenic. Assuming the radiative intensity to be given by the 
Planck function, the rate of production of electrons with energy in 
the range £ + d£ from a group of excited levels having energy in 
tl> ange X, X + uX below the ionization limit is 


n 

e 


I (2n m kT ) 
e e 7 


W 



-C/kT e 

dXd£ 

x + 5 


(C.6) 


The rate of production of electrons with energy in the range 
C, C + d£ from all the excited levels is 


e 2(Pnm kT ) 5/ 
x e e 7 


TT- 


Sh^H 


* -C/kT 

e 


/ dX 

J * ~ l 


'Therefore the avei^age energy of an electron created by photo- 
ionization from the excited states of the atom is 
“ -C/kT 


i 


In 


*M + ^ ' 


C 


J ^ 


l = 


° -C/kT f *M + ^ \ r 

b/ e in ! = — ! dc 

e V C J 


o 
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(C.7) 

Numerical integration of the integrals in equation (C.7) for the 

representative value X.7 kT = 4 gives 

iA e 

X = .64^ kT e (C.6) 

The average energy of electrons created in the atom-atom ionization 
process or by photoionization is of the order kT g or les as was stated 
in Chapter II. 
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APPENDIX D 

CALCUIATIONAL DETAILS 


As explained in Chapter II, the radiative terms are made arti- 
ficially small by a multiplying factor. Such a factor multiplies each 
of <4^, and <4^. Starting with the radiationless solution 

the factors for and may be immediately set to one since they 
have little effect on the solution behind the shock discontinuity. The 
factors for 14 and Q' must be small, however, since these terms 
control the radiative cooling and ionization rate. The procedure used 
was to keep the factor for small (usually zero) until the solution 
was found for which the factor for was one. In this way the long 

quasi-equilibrium region caused by the radiative ionization rate op- 
posing the collisional ionization rate near equilibrium need not be 
calculated each time, resulting in a considerable reduction in computer 
time. 

Successive steps in the solution are shown in Figure D.l where 
the degree of ionization in the radiative cooling region is shown. The 
degree of ionization is only slightly affected by the excited state 
radiation elsewhere. Curve number 1 is the solution obtained with the 
(o* u factor .1 and using as the previous solution that obtained with the 

OM 

14 „ factor zero. Curves 2 ana 5 show the solutions obtained when the 
oM 

Q u factor is increased to .3 and then to .5 each time using the previous 

oM 
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solution to estimate the radiative terms. (The larger value of the 
degree of ionization for curve 2 vbb caused by an error in the input 
values to the computer program, subsequently corrected.) The factor 
was then held at .5 and convergence was rapidly obtained resulting to 
curve 4. Using this solution and increasing the Q qM factor to .7 
results in curve 5* The Q u factor was then made one and tne average 
of the latest solution arvl the previous solution was used to generate 
the next solution. Convergence was rapid. The factor for was 

then raa-le one. The solution was little changed except for the delay 
of equilibrium which occurred. The result is shown as curve 6. 

It required about thirty minutes of IBM 70^4 computer time to 


reach t.he final solution. 


i 


